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! Abstract: 
> ■ 

in ; 

Tj- , For any arbitrary algebraic curve, we define an infinite sequence of invariants. We 

^ ! study their properties, in particular their variation under a variation of the curve, and 

^ . their modular properties. We also study their limits when the curve becomes singular. 

O ■ In addition we find that they can be used to define a formal series, which satisfies 

■ formally an Hirota equation, and we thus obtain a new way of constructing a tau 

^ I function attached to an algebraic curve. 

I These invariants are constructed in order to coincide with the topological expansion 

^ ! of a matrix formal integral, when the algebraic curve is chosen as the large N limit of 

>• . the matrix model's spectral curve. Surprisingly, we find that the same invariants also 

^ ■ give the topological expansion of other models, in particular the matrix model with an 

' external field, and the so-called double scaling limit of matrix models, i.e. the (p, q) 

minimal models of conformal field theory. 

As an example to illustrate the efficiency of our method, we apply it to the Kont- 
sevitch integral, and we give a new and extremely easy proof that Kontsevitch integral 
depends only on odd times, and that it is a KdV r-function. 
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1 Introduction 

Computing the topological expansion of various matrix integrals has been an interesting 
problem for more than 30 years. The reason for it, is that formal matrix integrals (cf 
[4T] for a definition of formal integrals), are known to be combinatorics generating 
functionals. Some formal matrix integrals enumerate maps, or colored maps of given 
topology [m [2ni EZ] , some count intersection numbers (the Kontsevitch integral and 
its generalizations [58]), and physicists have also tried to reach the limit of continuous 
maps through critical limits, and thus tried to recover Liouville's field theory [571 128] . 

Many methods have been invented to compute those formal matrix integrals, and 
the most successful is undoubtedly the "loop equations" method [53], El], which is 
in fact nothing but integration by parts, or Tutte's equations [73l or Schwinger- 
Dyson equations, or Ward identities, or Virasoro constraints, or W-algebra [26]. Until 
recently, those loop equations were solved only for the first few orders (mostly planar 
or torus), and case by case (for each matrix model). One of the most remarkable 
methods was obtained in [7]. Let us also mention that other methods were invented 
using orthogonal polynomials [67] (only in the case were the formal integral comes from 
an actual convergent integral), or topological string theory methods [651 [22l [11] using 
the so-called holomorphic anomaly equations. 

In 2004, a new method for computing the large expansion of matrix integrals 
was introduced in [33], and further developped in [iQl [17]. The starting point of that 
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method was not new, it was the same as in [7], it consists in solving the loop equations 
recursively in powers of the expansion parameter l/N"^, where is the size of the 
matrix. To leading order, loop equations become algebraic equations, and give rise to 
an algebraic curve S{x,y) = where S is some polynomial in 2 variables, which we 
call the "classical spectral curve" . 

The new feature which was introduced in [33j was to use contour integrals and 
functions on the curve rather than on the x-plane as in [7j. When written on the 
curve, the loop equations, together with the Cauchy residue formula and the Riemann 
bilinear identity, simplify enormously, and take a very universal structure which can be 
written entirely in terms of geometric properties of the curve. In other words, 
the solution of loop equations of many different matrix models, depends only on the 
properties of the spectral curve, and not on the matrix model which gives that curve. 
In particular, they can be written for any arbitrary algebraic curve, even for curves 
which don't come from matrix models. It is thus tempting to define "free energies" for 
any algebraic curve. This is what we do in this article. 

Therefore, in this article, for any arbitrary algebraic curve S{x, y) = 0, we define an 
infinite sequence of complex numbers F'^^\£), computed as residues of meromorphic 
forms on the curve. Out of these F^^\£)^s, we build a formal power series: 

oo 

InZjv(^) = -Y,N^'^'F^'\^) (1-1) 

9=0 

and we study its properties. 

We compute the variations of F^^^ under variations of the curve (variations of its 
complex structure, its moduli, and modular transformations). We show that the F^^^^s 
are invariant under some transformations of the curve, namely under transformations 
of the curve which preserve the symplectic form up to a sign ±dx A dy. 

We also show that Zn{S) satisfies bilinear Hirota equations, and thus Zi\f{S) is a 
formal r-function, and we construct the associated formal Baker-Hakiezer function [9] . 

We thus have a notion of a r function associated to an algebraic curve. Such 
notion has already been encountered in the litterature [I9j, and it is not clear whether 
our definition coincides with other existing definitions. What can be understood so 
far, is that we are defining a sort of quantum deformation of a classical r-function 
whose spectral curve is S. The classical r function being only the dispersionless limit 
In Zoo (i^) = —F^^\£), while our Zn{8) concerns the full system. 

Almost by definition, if £ is the algebraic curve coming from the large N limit of 
the loop equations of a matrix model, then Zi^[£) is the matrix integral. 
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What is more interesting is to see what is Zn{S) for curves not coming from the 
large N hmit of the loop equations of a matrix model. 
We study in details a few examples. 

- The double scaling limit of a matrix model. It has been well known since [5T| [26]. 
that if we fine tune the parameters of a matrix model so that the algebraic curve S 
develops a singularity, the free energies become singular and the most singular part of 
the free energies form the KP-hierarchy r-function (KdV hierarchy for the 1-matrix 
model). We show, by looking at the double scaling limit of matrix models, that the 
KP r-function (resp. KdV r-function), coincides with our definition for the classical 
limit of the {p,q) systems (resp. (p, 2)). 

- It has been well known since the works of Kontsevitch [58j, that the KdV r- 
function can be represented by another matrix integral called Kontsevitch integral. 
Kontsevitch introduced that integral as a counting function for intersection numbers, 
and proved that it is a KdV r-function. One of the key features is that it depends on 
the eigenvalues of a diagonal matrix A, only through the quantities tk = tr for 
odd k (cf lini EH]). Another important known property is that, if tfc = for k > p, it 
coincides with the {p, 2) r-function found from the double scaling limit of the 1-matrix 
model, i.e. the {p, 2) conformal minimal model. 

Here, we prove that the Kontsevitch matrix integral coincides with our Zn{S) when 
S is the large N limit of the Schwinger-Dyson equation of the Kontsevitch integral. 
The remarkable fact, is that for our Z]\f{S), the above properties (i.e. the fact that it 
depends only on odd t^s and the fact that it gives the (p, 2) r-function if = for 
k > p) are trivial. We thus provide a new proof of those properties, and maybe a new 
interpretation. 

2 Main results of this article 

In this section we just sketch briefly the contents of the main body of the article. 
2.1 Definitions 

Given a polynomial of two variables S{x,y), we construct an infinite sequence of mul- 
tilinear meromorphic forms over the curve of equation S{x,y) = 0, which we call: 

Wt\pi,p2,...,pk) , k,geN. (2-1) 

In particular Wq^^ = —F^^^ are complex numbers F^^\S). 

The F(f)'s and the VT^^^'s are defined in terms of residues near the branchpoints of 
the curve only, i.e. they depend only on the local behavior of the curve near its 
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branch points. 

Then we show some properties: 

• the Wl^^^'s are symmetric in their k variables; 

• there is a "loop insertion operator" which increases k ^ k + 1: 

DBip,+„.)W^^\pi,P2, ■■■,Pk)^ W^fe?i(Pi,P2, ■ ■ ■,Pk,Pk+i); (2-2) 

• there is an inverse operator which contracts k ^ k — 1: 

Res . ^Pk)wi'\p,,p2,...,Pk)^{2g + k-3)wli'\{p,,p2,...,Pk-i). 

Pfc-^branch points 

(2-3) 

The F(5)'sand theW^^^'s are defined in a way which mimics the solution of matrix 
models loop equations, and almost by definition, they coincide with matrix model's 
free energy and correlation functions when the polynomial £ is chosen as the classical 
large N limit of the matrix model's spectral curve: 

ln( f dM exp -N Tr V{M)) = -f^ ^^'^^ ^^^^ {^imm) (2-4) 



and 

^ ■ ■ ■ ^ ■ ■ ■ 

(2-5) 

The same construction works also for the 2-matrix model and the matrix model in 
an external field: 

Hum = F^^\^2Mm) 

^e2.field = ^^'n^ext.field). (2-6) 

in particular it works for the Kontsevitch integral 

-^Kontscvitch ~ (^Kontsevitch) (2-7) 

where the LHS is the topological expansion of the corresponding matrix integral, and 
the RHS is the functional defined in this article, applied to the curve S{x, y) = coming 
from the large A*" limit of the Schwinger-Dyson equations of the corresponding matrix 
model. 

Let us emphasize that not every curve S is the large N limit of a matrix model's 
spectral curve, and thus our functional F^^\S) is defined beyond matrix models, and 
is really an algebro-geo metric object. It has many remarkable properties, and we list 
below some of the most important ones: 
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2.2 Remarkable properties 
Theorem 14.81 Diagrammatic representation: 

/ \ 



wl^_l-^{p,Pi,...,Pk) = Yl w{G)=w 



Gegi%{p,pi,...,p^,) 



(2-8) 



where Q'f_li{p, pi, . . . , pk) is a set of trivalent graphs (built on trees), and w is a Feynman- 
like weight function associating values to edges and integrals (residues in fact) to ver- 
tices of the graph. 

This theorem is important because it makes all formulae particularly easy to re- 
member, and many theorems below can be proved in a diagrammatic way. 

Theorem 18.11 Singular limits: if the curve becomes singular, the functional F^^-* 
commutes with the singular limit, i.e.: 

lim ) {S) = ) (lim S). (2-9) 
Theorem 19.21 Integrability: the formal series: 

oo 

In Zn{S) = -J2n^~^'F^'\^) (2-10) 

9=0 

obeys Hirota's bilinear equations, and thus is a r function. 

Theorem 15.31 Homogeneity : F^^\£) is homogeneous of degree 2 — 2g in the moduli 
of the curve: 

gpig) 

(2-2^)F(^) = ^4-— . (2-11) 

Theorem 15.11 Deformations : if the curve £ is deformed into £ + 5£, the differential 
ydx is deformed into ydx ydx + 6{ydx) where S{ydx) is a meromorphic one- form 
which we denote 6{ydx) = — fi, and which can be written as: = J^^ W2^^A for some 
appropriate contour dfl and some appropriate function A. Then we have: 

SWt'\puP2, ...,Pk)= [ A{Pk+i) W^ifi(Pi,P2, . . .,Pk,Pk+i). (2-12) 

Jan 

Theorem 17.11 Symplectic invariance : F^^\£) is unchanged under the following 
changes of curve £{x, y): 

y ^ y + R{x) , _R(x)=rational fraction of x, 

y ^ cy , X ^ c~^x , c=complex number, , . 

y ^ -y , X ^ X, 

y ^ X , X ^y. 
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all those transformations conserve the symplectic form dx A dy up to the sign. 

Theorem 16.21 Modular transformations : the modular dependence of F^^\8) 
is only in the Bergmann kernel (defined in section 13.1.51) . and thus the modular 
transformations of F^^\£) are derived from those of the Bergmann kernel. Under 
a modular transformation, the Bergmann kernel is changed into B{p,q) B{p,q) + 
2iTT du{p) Kdu{q) , and thus we introduce a new kernel for any arbitrary symmetric ma- 
trix K,: 

B^{p,q) B{p,q) + 2i7rdu{p)Kdu{q) (2-14) 
We thus define some F!f\S), and we compute: 

We also remark that when k = (r — r)~^, F!f\8) is modular invariant. 
2.3 Some applications, Kontsevitch's integral 

{p, q) minimal models, KP and KdV Hierarchies: it is well known [26l |20] that 
some rational singular limits of matrix models correspond to {p, q) minimal models, 
and theorem 18.11 implies that: 

fS.^ = F^'\£m) (2-16) 

and it is well known that {p, q) minimal models are some reductions of KdV hierarchy 
for q = 2 and KP hierarchy for general {p, q). 

Notice that the x ^ y symmetry of theorem 17. II (i.e. Eq. fl2-13p ) implies the famous 
{p,q) ^ {q,p) duality [26l [281 [56] . 

Kontsevitch integral's properties: Kontsevitch's integral is defined as: 

■Z^Kontsevitch(A) = dM e ^'^^ » , lu .^Kontsevitch = — N'^ -^Kontsevitch 

9=0 

(2-17) 

and we define the Kontsevitch's times: 

4 = ^TrA-^ (2-18) 

It is straightforward to write the Schwinger-Dyson equations and find the classical 
spectral curve: 

x{z) = z + j^Ti -T—T 



P.r _ ) ■^K'^l 2iV A^-A ('9-1 Q^l 

'-'Kontsevitch l/^ 2ia. ■ \ ) 

' y[z) = + ti 
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According to theorem I10.3[ we have: 

-^Kolitsevitch = -^^^n^Kontsevitch)- (2-20) 

Using the x ^ y invariance of eq J2-13[ we see that the only branch point in y is located 
at 2; = 0, and since the F^^^^s depend only on the local behavior near the branchpoint, 
we may perform a Taylor expansion of x{z) near z = 0: 

t^Kontsevitchltl, 12, ^^^j ^ ^2 _^ • l^-^J^J 

From the symplectic invariance theorem 17.11 (i.e. eq J2-13l) . we may add to x any 
rational function of y i.e. of thus we may substract to x its even part, and thus the 
following curves are related by symplectic invariance: 

^Kontsevitch (^1, ^2-, • • •) ~ ^Kontsevitch (if^l, 0, ^3, 0, ^5, . . .). (2-22) 

We thus have a very easy proof that -^Kontsevitch depends only on odd times. 

Moreover, if = for A; > p + 2, we have: 

^Kontsevitch (^1, ^2, • • • , ^p+2, 0, . . .) = | _ ^2 .j.^^^'^"" ^'^'^ (2-23) 

which is exactly the curve of the (p, 2) minimal model, i.e. it satisfies KdV hierarchy. 
We thus have a very easy proof that ^'xontscvitch is a KdV tau-function. 

Those are old and classical results about the Kontsevitch integral, and we just 
propose a new proof, in order to illustrate the power of the tools we introduce. 

3 Algebraic curves, reminder and notations 

We begin by recalling some elements of algebraic geometry, which are used to fix the 
notations. We refer the reader to [46j or [47j for further details about algebro-geometric 
concepts. 
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Summary of notations 



S{x,y) = — > classical spectral curve. 

di + 1 = deg^ S — > x-degree of the polynomial S. 

d2 + I = degy £ — > y-degree of the polynomial £ (number of sheets). 

a — {oj} — > set of branch points dx{ai) ~ 0. 

a = {ccj} — > poles of ydx. 

^ — > genus of the curve. 

Ai n Bj — Sij — > cannonical basis of non-contractible cycles. 

dui — > cannonical holomorphic forms dui = 6ij. 

Tij = §„ dUi — > Riemann's matrix of periods. 

—j 

Uiip) = dui — > Abel map. 

A — A — k{B — tA) — > K-modified ^-cycles 

B — B — tA — > K-modified B-cycles, Ai n Bj — Sij 

dSq^^q^{p) — > 3rd kind differential with simple poles qi and q2, 

such that ReSgj dSq^^q^ — 1 — — Res q^dSq^^q^ and 

B{p, q) — )• Bergmann kernel, i.e. 2nd kind differential with 

double pole at p — q, no residue and vanishing A- 
cycle integrals. 

^ _ n+rm — ^ regular odd charactcrisic, i.e. Yli=i =odd. 

dhz — ^^r^ d,Ui — > Holomorphic form with only double zeroes. 



i)=0. 



K(P: q) — I — Prime form mdependent of with a simple zero at 

p = q. 

= ydx — > some antiderivative of ydx defined on the universal 

covering. 

p, x{p) — x{p) — > if p is near a branchpoints a, then p ^ p is the 

unique other point near a such that x(p) — x{p). 

p^{p)i x{p^) = x{p) — > The p*'s, i = 0, . . . ,d2, are the pre-images of x{p) 

on the curve. By convention p^{p) = p. 

Dft = Sq + ti {k,6qt Kj^) — »• Covariant variation wrt Q — 6{ydx). 



Consider an (embedded) algebraic curve given by its equation: 

£{x,y) = (3-1) 

where £ is an almost arbitrary polynomial of two variables. This is equivalent to 
considering a compact Riemann surface E and 2 meromorphic functions x and y, such 
that 

VpeE , £(x{p),y(p))^0. (3-2) 

We only require that £{x,y) is not factorizable, and that all branchpoints (zeroes 
of dx) are simple, i.e. near a branchpoint a^, y behaves like a square root ■\/ x — x{ai). 
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x(q) x(p) X 



Figure 1: Example of an algebraic curve with two x-brancli points Oi and 02 and a 
three sheeted structure {x has three preimages). One can see that the map p —>■ pis not 
globally defined, for instance when g — > p, we have q p^'^\ The notion of conjugated 
point depends on the branch point. 

3.1 Some properties of algebraic curves 

3.1.1 Sheets 

For each complex x, there exist d2 + 1 = degyE solutions for y of S{x,y) = 0. This 
means that there are exactly d2 + l points on the Riemann surface S for which x{p) = x: 
E has a sheet structure with 6/2 + 1 x-sheets. We call them: 

x{p) = X <-> P = p\x) , i = 0,...,d2- (3-3) 

3.1.2 Branch points and conjugated points 

Let Qi, i = 1, . . . ,n, SL = {ai, . . . , a„} be the set of branch-points, solutions of dx = 0: 

Va G a , dx{a) = 0. (3-5) 

Since we assume that the branch-points are simple zeros of dx, we have the following 
property: if p is in the vicinity of a branch-point Oj, there is a unique point p ^ p, such 
that x{p) = x{p), which is also in the vicinity of Oj. p depends on i, and in general, p 
is not globaly defined (see fig. [Hfor an example). 

Notice that p is one of the defined in the previous section. 

3.1.3 Genus, cycles, Abel map 

If the curve has genus ^ , there are 2^ homologicaly independent non-trivial cycles, 
and we may choose a (not unique) cannonical basis: 

Ai n = 6ij , n 4. = , Bif] Bj = 0. (3-6) 
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(3-11) 



Figure 2: Example of canonical cycles and the corresponding fundamental domain in 
the case of the torus (genus ^ = 1)- 



The simply connected domain obtained by removing all A and S-cycles from the curve 
is called the "fundamental domain" (see fig. ([2]) for the example of the torus). 

On a genus ^ curve, there are ^ linearly independent holomorphic forms 
dui, . . . , du^ , which we choose normalized on the ^-cycles: 



-4. 



dui = 6ij. (3-7) 

The Riemann matrix of period is defined by the S-cycles 

Tij = (t dui. (3-8) 



They have the property that 

Tij = Tji , lmr>0. (3-9) 

Given a base point pq on the curve (we assume it is not on any A or ^-cycle), we 
define the Abel map 

Ui{p) = / dui (3-10) 
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where the integration path is in the fundamental domain. The ^ -dimensional vector 
u{p) = {ui{p), ■ ■ ■ (p)) maps the curve into its Jacobian. 

3.1.4 Theta-functions and prime forms 

We say that z G is a characteristic if there exist two vectors with integer coeffi- 
cients a E and b E such that: 

d + T.b 

z = . (3-12) 

z is called an odd characteristic if 

^^aj6j = odd. (3-13) 

i=l 
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Given a characteristic z = ^^j-^, and given a symmetric matrix Tij = r,j such that Imr 
is positive definite, and a vector v E , we define the 9^ function: 

^ gj7r(n-6/2)'r(n-6/2) g2j7r(t.+a/2)t.(n+6/2)_ (3-14) 

Hez^ 

If z is an odd cliaracteristic, 6z is an odd function of v, and in particular 6^(0, r) = 
and we define the following holomorphic form: 



do (v^ 

dh,{p) = ^dui{p) . 



. , dvi 
1=1 



(3-15) 

't;=0 



All its ^ — 1 zeroes are double zeroes, so that it makes sense to consider its square 
root defined on the fundamental domain. The prime form is: 

E{ ) ^^(^(^') ~ ^(^)) ^2 j^g^ 

It is independent of z, and it vanishes only if p = g (and with a simple zero), and has 
no pole. 

3.1.5 Bergmann kernel 

We define a bilinear meromorphic form called the "Bergmann kernel" (TOlHH]: 

B_{p, q) = Bergmann kernel (3-17) 

as the unique 1-form in p, which has a double pole with no residue aX p = q and no 
other pole, and which is normalized such that 

mv.q) //l^^'^l^l^ + fi^^t^ ' I B = Q (3-18) 

[z[p) - z[q)y J^^ 

where z is any local coordinate on the curve in the vicinity of q. The Bergmann kernel 
depends only on the complex structure of the curve, and not on the details of £. For 
instance if £ has genus zero, B_ is the Bergmann kernel of the projective complex plane 
(the Riemann sphere), and if the curve has genus 1, 5 is related to the Weierstrass 
function. 

Properties: 

mp,q) = B{q,p) , / B{p,q) = 2mdUi{p). (3-19) 

For any odd characteristic z^ we have: 

B{p, q) = dpdgln {9,{u{p) - u(g))) (3-20) 
If f{p) is any meromorphic function, its differential is given by 

df{p)= Res Bip,q)f{q). (3-21) 



14 



3.1.6 3rd type differentials 

Given two points qi and q2 on the curve , we define the 1-form dS_q-^^q^ by: 

dSquJp)- r B{p,q) (3-22) 

where the integration path is in the fundamental domain. 

—91,92 unique meromorphic form with only simple poles at gi and q2i such 

that: 

Res dS^^^^^ = 1 = - Res dS^^^^^ , £ dS^^^^^ = 0. (3-23) 

3.1.7 Modified set of cycles: 

In order to easily deal with modular properties of the objects we are going to introduce, 
it is convenient to define some modified cycles and kernels with an arbitrary symmetric 
matrix k. When k = 0, all those quantities reduce to the unmodified ones. The 
modular transformations of the modified objects, merely amount to a change of k. 

We thus choose an arbitrary ^ x ^ symmetric matrix k with complex coeffi- 
cients, and we define another set of cycles: 

j I 

B, = B,-Y,r^,A, (3-24) 

3 

They satisfy: 

Ai n Bj = , A n A- = , n B^- = (3-25) 

and we straightforwardly have 

duj = 6ij , ® duj = 0. (3-26) 

3.1.8 Modified Bergmann kernel 

We also define the modified Bergmann kernel, normalized on A instead of A: 

B{p, q) = B{p, q) + 2m ^ dui{p) Hij duj{q). (3-27) 

i,3 

It is such that 

B{p,q)^B{q,p) , (f B ^ , (f B{p,q) ^ 2m du.ip) (3-28) 

JAi J qeBi 
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and if f{p) is any meromorphic function, its differential is given by: 

df{p)= Res B{p,q)f{q). (3-29) 

• For K = we have B_ = B. 

• For K = (r — t)~^, B is the Schiffer kernel [TOl US], and it is modular invariant. 

3.1.9 Modified prime form 

Similarly we define a modified prime form: 

E{p, q) = E{p, q) e2^-«*(p) (3_30) 
It vanishes only ii p = q (with a simple zero), and has no pole. 

3.1.10 Modified 3rd type differentials 

In the same fashion, we define the modified 3rd type differentials dSq^^q.^ by: 

dSq,,qM= r B{p,q) (3-31) 
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where the integration path is in the fundamental domain. 

d'Sqi,q2 is the unique meromorphic form with only simple poles at qi and q2, such 
that: 



Res dSq^ 52 = 1 = — Res dSq^ q^ , (f dSq^ q2 = 0. (3-32) 
Properties: 

d^qi,q2 ~ ~d^q2,qi (3-33) 

dSq,,qM = dpln (^ g^l^^l + 2ivr dui{p)Kij{uj{qi) - Uj{q2)) (3-34) 

f dSq^^q2 = 2m{ui{qi) - Ui{q2)) (3-35) 

dq, {dSq,,q,{p)) = B{qi,p) (3-36) 

I dSq^^q^ = / dSp^^p^ (3-37) 
Jpi J qi 

Cauchy residue formula: for any meromorphic function f{p) we have: 

f{p) = -ResdSq,^qMf{qi)- (3-38) 

91 
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3.1.11 Bergmann tau function 

The Bergmann r-function tbx was introduced and studied in [6T | [621 [37], it is such 
that _ 

?^ = R« ^ (3-39) 

ox[ai) p^ai dx[p) 

It is well defined because the Ranch variational formula [68j implies that the RHS is a 
closed form. Notice that tbx is defined only up to a multiplicative constant which will 
play no role in all the sequel. 

3.2 Examples: genus and 1 

• Genus 

If the curve S has a genus ^ = 0, it is conformaly equivalent to the Riemann sphere, 
i.e. the complex plane with a point at oo, and there exists a rationnal parametrization 
of the curve. It means that there exists two rational functions X{p) and Y{p) such 
that: 

£{x, y)=0 ^ 3peC, x = X{p) , y = Y{p) (3-40) 
In this case, the Bergmann kernel is the Bergmann kernel of the Riemann sphere: 

B{p, q) = B{p, q) = = d, d. In (p-q). (3-41) 

The prime form is: 

E{p,q) = E{p,q) = ^^. (3-42) 
\/dpdq 

• Genus 1 

If the curve has genus ^ =1, then it can be parametrized on a rhombus corresond- 
ing to the fundamental domain of a torus (see fig. [2]). It means that there exists two 
elliptical functions X{p) and Y{p) such that (see [75] for elliptical functions): 

£{x, y) = ^ 3p ^ C , X = X{p) , y = Y{p) 
X{p+l) = X{p + T)=X{p) , Y{p + 1) =Y{p + t) =Y{p) (3-43) 

Then, the Bergmann kernel is the corresponding Weierstrass function [75] : 

B{p, q) = {p{p -q,r) + ^) dpdq. (3-44) 

The prime form is: 

When K = 2ii^r ' modified Bergmann kernel is the Schifii^er kernel, and if ^ = 1 
it is the Weierstrass function: 

B{p, q) = p{p - g, r) dpdq. (3-46) 
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3.3 Riemann bilinear identity 

If uji and UJ2 are two meromorphic forms on the curve. Let po be an arbitrary base 
point, we consider the function $i defined on the fundamental domain by 

= Tcui. (3-47) 

We have 

Res $i(p)cj2(p) = ^— ^ / ui (f UJ2 — (f uji (f UJ2- (3-48) 

p^all poles 2z7r ^ J^, J B, J B^ J 

Note that this identity holds also for the modified cycles with any k: 



Res $1(^)6^2 (p) = ^— ^ / ui (f UJ2 — (f uji (f U2- (3-49) 

p^all poles 2«7r ^ Ja J J B, J Ai 

In particular with uJi{p) = B{p, q), we have: 



Res dSp,po{q) a;(p) = - V rfMi(g) * 

p—> all poles ^— ' / < 



u (3-50) 



and: 



;(g) = Res (iS'p.po(g) c^(j9) + ^^^(g) cu. (3-51) 

p^polesofw ' ' J A 



p^poles oi oj 

1=1 

3.4 Moduli of the curve 

The curve £{x, y) = is parameterized by: 

• a genus ^ compact Riemann surface S, with periods Xjj. 

• punctures at the poles of x and y, whose moduli are given by the negative 
coefficients of the Laurent series of ydx near the poles. 

• the ^-cycle integrals of ydx, called filling fractions. 

3.4.1 Filling fractions 

We define: 

e^ = ^ (f ydx (3-52) 
which are called "filling fractions" by analogy with matrix models (|35j). 
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3.4.2 Moduli of the poles 

Consider a pole a of ydx^ define the "temperatures" : 

tot — Res ydx. (3-53) 

a 

Notice that: 



5^ta = 0. (3-54) 



Then consider the 3 cases: 

• Either ck is a pole of x of degree da-, then we define the local parameter near a as: 

Zcx{v)=Av)'^\ (3-55) 

• or q; is not a pole of x neither a branchpoint (thus it is a pole of y), then we define 
the local parameter near a as: 

^a(p) = , , y (3-56) 
x\p) — xya) 

• or a is not a pole of a;, and it is a branchpoint (thus it is a pole of y), then we 
define the local parameter near a as: 

= , ^] , (3-57) 

In all cases, in the vicinity of a, we define the "potential" 

K(p) = Res y{q)dx{q) In U - (3-58) 
which is a polynomial in Za.{p): 

deg Va 

K(p)- W^a(p)- (3-59) 

fc=l 

The ^Q,*; are the moduli of the pole a. 

We have the following properties: 

dVM = Res y{q)dx{q) (3-60) 

Res dV^ = (3-61) 



a 



and 

,(p).xW^^^dK„(p)-t„i||;l + 0(|^). (3-62) 
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We have from Eq. (13-511) : 

y{p)dx{p) = - Res B{p,q)Va{q) + y^tadSa,o{p) + 2iTT eidui{p) . (3-63) 



If we introduce 



Bo.Av) = - Res B{p, q) z^{q)\ (3-64) 



we can turn this expression to 

ydx = ^ ta,kBa,k + ^ tadSa,o + 2ivr ^ eidui (p) (3-65) 

a,k a i 

in order to exhibit the moduh of the curve. 

4 Definition of Correlation functions and free ener- 
gies 

In all this section, the curve S{x, y) = and a symmetric matrix k are given and fixed. 
The unfamiliar reader may choose n = since most usual applications (matrix models) 
correspond to that case. 

4.1 Notations 

Consider an arbitrary point p G E, and a point g of E which is in the vicinity of a 
branch point aj (so that q is well defined). We define: 

Definition 4.1 Diagrammatic rules: 

vertex : uj{q) = {y{q) - y{q))dx{q) (4-1) 

line — propagator : B{p, q) (4-2) 

1 

arrow — propagator : dEq{p) = - i B{^,p) (4-3) 

2 Jq 

where the integration path is a path which lies entirely in a vincinity of ai ( thus it is 
uniquely defined^. 

root : ^{q) = ydx (4-4) 

J o 

where o is an arbitrary base point on the curve, i.e. $ is an arbitrary antiderivative of 
ydx, i.e. d^ = ydx. 



■^This definition is the opposite of the notation used in [40l [19] since the integral goes from q to q 
instead of going from q to q. 
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The reason why we call these objects diagramatic rules and vertices, propagator or 
root, is explained in section H75l below. 

Notice that dE depends on i, i.e. on which branchpoint we are considering, but we 
omit to mention the index i in order to make the notations easier to read. In all what 
follows it is always clear which i is being considered. 



Notation for subset of indices: 

Given a set of points of the curve {pi,P2, ■ ■ ■ ,Pn}, if K 
subset of {1,2,..., n}, we denote: 

PK = {Pil,Pi2, ■ ■ • ,PiJ 



{ii,i2, ■■■,ik} is any 



(4-5) 



4.2 Correlation functions and free energies 
4.2.1 Correlation functions 

The fc-point correlation functions to order g, wlf \ are meromorphic multilinear 
forms, defined by the following recursive triangular system: 

Definition 4.2 Correlation functions 



W, 



if ^ < 



= 



wi'\pi,p2)=B{p,,p2) 



and define recursively (remember that pk is a k-uplet of points cf eq \4-5\ ): 



H^ii(p,PK)= Res, 



il|j7ii(9, Pk/j) + W^'i2^\q, q, Pk 



(4-6) 
(4-7) 

(4-8) 



(4-9) 



This system is triangular because all terms in the RHS have lower 2g + k than in the 
LHS and given Wi^^ and W2^\ it has a unique solution. 

Notice that wlf_^^{p,pi, . . . ,pk) is a multilinear meromorphic form in each of its 
arguments, it is clearly symmetric in the last /c-ones, and we prove below (theorem 14.61) 
that it is in fact symmetric in all its arguments. 

More properties of wjf^-^ are studied below in section 14.41 
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4.2.2 Free energies 

We define tlie free energies which are complex numbers: 

Definition 4.3 Free energies. 
For g > I 




and 




where 



y [at 



dy{ai 



Zi{p) = \/x{p) - x{ai) 



dzi{ai) 

and Tbx is the Bergmann r-function defined in Eq. h3-39^ . 
F^^'^ is defined in the next section. 

4.2.3 Leading order free energy F^^K 
Let us define F^^^ as follows : 



(4-10) 



(4-11) 
(4-12) 



= ^ Res Vaydx + ^ t^fJ'a - f f ydx 



(4-13) 
where /i^ is given by 

r dz 
fi^= / iydx-dV^ + t^—) + V^{o)-tJn{z^{o)) (4-14) 

J a 

Notice that /Iq depends on some base point o, but the sum ^^^ta/ia does not. 
4.2.4 Special free energies and correlation functions 

All the quantities defined so far, were defined with the K-modified cycles and modified 
Bergmann kernel. Let us also define them for k = (for instance F^^^ obviously needs 
another definition). 

Therefore we also define the unmodified quantities corresponding to /t = 0, as: 



^k,g, Wl'\p^,...,pk)■.= Wt'>{pu...,Pk) 
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rig) I 



K=0 



(4-15) 



for g > 1, := -i— V Res ^q)wi'\q) = , (4-16) 

F^'^ = -\MrB.) - ^In {\{y'{a,)^ , (4-17) 

and = - Res Vaydx + - ta^ia — / ydx (p ydx . (4-18) 

2^ a 2 ^ 4z7r ^ /_4. J IS 

Remark 4.1 The special functions, except F_^^^ and F_^'^\ are obtained by changing B and 
dS by S and dS_ in the diagrammatic rules defined in section I4.5[ 



4.3 Tau function 

Definition 4.4 We define the tau-function as the formal power series in N^"^: 



9=0 

(4-19) 

We sliow in section [9] that Z]\r{S) is indeed a tau-function because it obeys Hirota 
equations, order by order in A^~^. 

4.4 Properties of correlation functions 

The loop functions defined in definition 14.21 satisfy the following theorems, whose 
proofs can be found in Appendix [Al 

Theorem 4.1 The correlation function W^^^ is worth: 

w(o)^ x o B{q,p)B{q,pi)B{q,p2) 

W^, \p,p,,p,) = Res ^-^-^^ (4-20) 

In particular, ^^3'''' is symmetric in its 3 variables. 

Theorem 4.2 For {k,g) 7^ (1,0), the loop function wl:^_^^{p,pi, . . . ,pk) has poles (in 
any of its variables p,pi, . . . ,pk) only at the branch points. 

Theorem 4.3 For every A cycle we have: 

Vz = l,...,^ / W^ifi(p,pi,...,pfc) = 0, (4-21) 

Vz = l,...,^ ,Vm = l,...,fc / W^ifi(p,pi,...,pfc)=0. (4-22) 
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Theorem 4.4 For every k and g, we have: 

p\pi,...,Pk) _ 

dx{p^) "'^ ix{p) — x{pi)) 

and if k > 1: 



^ W^'i^{p\pi,...,Pk) _ dxjpi) 

2^ (^(^■._^(^^2 ^^-^"^l 



sr^ wifl^{pi,p\P2,...,Pk) _rr dxjpi) 
^ dx{p^) ~ {x{p) - x(pi))2 

where we recall that p^ are all the points such that x{p^) = x{p) (see section l3J\) . 

Theorem 4.5 For {k,g) ^ (0,1), 

Pt\xip), Pk) = [ - Up')dx{p)wt^,ip\ pk) + W^,%'\p\p\ pk) 



m=0 JCK 

(4 - 25) 

is a rational function of x{p), with no poles at branch- points. 
Theorem 4.6 W^'^^ is a symmetric function of its k variables. 



Corollary 4.1 



Vi, Resiyif,(p,pi,...,pfc) = 0, (4-26) 



Vz, Res x{p)wl:'^,{p,pu . . . ,pfc) = 0, (4-27) 

V Res |/(p)W^ifi(p,Pi, ...,Pk) = 0, (4-28) 

i 

V Res xip)y{p)wl^'^,ip,p,, . . . ,Pfc) = 0. (4-29) 



Theorem 4.7 For k > 1 we have: 

Res $(Pfc+i)Vrifi(PK,Pfc+i) = (2(7 + A; - 2)W^i^)(pK) + Sg,o5kMPi)dxiPi). 

(4-30) 

Notice that for k = and (7 > 2, it holds by definition if we define W^^^ = —F^3\ 
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4.5 Diagrammatic representation 

The recursive definitions of wjf^ and F^s) 

can be represented graphically. 

We represent the multihnear form wl^\pi, . . . ,pk) as a blob-hke "surface" with g 
holes and k legs (or punctures) labeled with the variables pi, . . . and F^^'^ with 
legs and g holes. 




(4-31) 

We represent the Bergmann kernel B{p,q) (which is also W2^\ i.e. a blob with 2 
legs and no hole) as a straight non-oriented line between p and q 

B{p,q) := p q . (4-32) 

We represent as a straight arrowed line with the arrow from p towards q, and 

with a tri-valent vertex whose legs are q and q 




(4-33) 



Graphs 

Definition 4.5 For any k > and g > such that k + 2g > 3, we define: 

Let ^j^+i(p,Pi, . . . ,Pk) be the set of connected trivalent graphs defined as follows: 

1. there are 2g + k — 1 tri-valent vertices called vertices. 

2. there is one 1-valent vertex labelled by p, called the root. 

3. there are k 1-valent vertices labelled with Pi, ■ ■ ■ ,Pk called the leaves. 

4. There are 3g -\- 2k — 1 edges. 

5. Edges can be arrowed or non-arrowed. There are k -\- g non-arrowed edges and 
2g -\- k — 1 arrowed edges. 

6. The edge starting at p has an arrow leaving from the root p. 
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7. The k edges ending at the leaves pi, . . . , are non-arrowed. 

8. The arrowed edges form a "spanning plana)^ binary skeleto^ tree" with root p. 
The arrows are oriented from root towards leaves. In particular, this induces a 
partial ordering of all vertices. 

9. There are k non-arrowed edges going from a vertex to a leaf, and g non arrowed 
edges joining two inner vertices. Two inner vertices can be connected by a non 
arrowed edge only if one is the parent of the other along the tree. 

10. If an arrowed edge and a non-arrowed inner edge come out of a vertex, then the 
arrowed edge is the left child. This rule only applies when the non- arrowed edge 
links this vertex to one of its descendants (not one of its parents). 

We have the following useful lemma: 
Lemma 4.1 There is a 1 to 3g -\- 2k — 1 map from ^^^^^(p, pk) to ^^+2(P' PK,Pfc+i)- 
proof: 

If G is a graph in ^f.+2(P)Pi) • • • iPkiPk+i)-, remove the non-arrowed edge attached 
to the leaf Pk+i and remove the corresponding vertex, and merge the incoming and the 
other outgoing edges of that vertex. You clearly get a graph G' G Pi, • • • iPk)- 

It is clear that the same graph is obtained 35^ + 2^ — 1 times (the number of edges of 
G' . And it is clear that from any G' G ^^.^^(p, Pi, • • • ,Pk), you can obtain 3(7 + 2 A; — 1 
graphs G G Q^^liiPiP^i ■ ■ ■ iPk.Pk+i) by adding a new vertex on any edge, and linking 
this new vertex to the leaf Pk+i- □ 

Example of Qf\p) 

As an example, let us build step by step all the graphs of Q{ (p), i.e. g = 2 and 
k = 0. 

We first draw all planar binary skeleton trees with one root p and 2g -\- k — 1 = 3 
arrowed edges: 

. (4-34) 

Then, we draw g + k = 2 non-arrowed edges in all possible ways such that every vertex 
is trivalent, also satisfying rule 9) of definition 14.51 There is only one possibility for the 




''it goes through all vertices. 

^planar tree means that the left child and right child are not equivalent. The right child is marked 
by a black disk on the outgoing edge. 

®a binary skeleton tree is a binary tree from which we have removed the leaves, i.e. a tree with 
vertices of valence 1, 2 or 3. 
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first graph and two for tlie second one: 





It just remains to specify the left and right chi 
possibilities in accordance with rule 10) of def J4.5l ard 



(4-35) 



dren for each vertex. The only 






(4-36) 

In order to simplify the drawing, we can draw a black dot to specify the right child. 
This way one gets only planar graphs. 






(4-37) 

Remark that without the prescriptions 9) and 10), one would get 13 different graphs 
whereas we only have 5. 

Weight of a graph 

Consider a graph G G ^^^^(p, Pi, • • • ,Pk)- Then, to each vertex i = 1, . . . ,2g + k — 1 oi 
G, we associate a label qi, and we associate qi to the beginning of the left child edge, 
and gj to the right child edge. Thus, each edge (arrowed or not), hnks two labels which 
are points on the Riemann surface S. 



To an arrowed edge going from q' towards q, we associate a factor 



{y{<i)-y{<i))d.x{q) ■ 

To a non arrowed edge going between q' and q we associate a factor -B(g', q). 



Note that the graphs are not necessarily planar. 
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Following the arrows backwards (i.e. from leaves to root), for each vertex g, we 
take a residue at g — a, i.e. we sum over all branchpoints. 



After taking all the residues, we get the weight of the graph: 

w{G) 

which is a multilinear form in p,pi, . . . ,pk. 

Similarly, we define weights of linear combinations of graphs by: 

w{aGi + (3G2) = aw{Gi) + (3w{G2) 

and for a disconnected graph, i.e. a product of two graphs: 

w{GiG2) = w{Gi)w{G2). 



(4-38) 



(4-39) 



(4-40) 



Theorem 4.8 We have: 



( 



Gegi%{p,pi,...,Pk) 



G 



(4-41) 



proof: 

This is precisely what the recursion equations 14-91 of def J4.2l are doing. Indeed, one 
can represent them diagrammatically by 





(4-42) 



□ 



Such graphical notations are very convenient, and are a good support for intuition 
and even help proving some relationships. It was immediately noticed after [33J that 
those diagrams look very much like Feynman graphs, and there was a hope that they 
could be the FejTiman's graphs for the Kodaira-Spencer theory. But they ARE NOT 
Feynman graphs, because Feynman graphs can't have non-local restrictions like the 
fact that non oriented lines can join only a vertex and one of its descendent. 

Those graphs are merely a notation for the recursive definition 14. 2[ 
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Lemma 4.2 Symmetry factor: The weight of two graphs differing by the exchange 
of the right and left children of a vertex are the same. Indeed, the distinction between 
right and left child is just a way of encoding symmetry factors. 

proof: 

This property follows directly from theorem 14.41 and the definition Eq. f l4-9p . Con- 
sider one term contributing to the first part of RHS of Eq. (14-91) : 



= - Res ^W;-Uq, Pj)Wt-iJ;Uq, Pk/j) (4-43) 
= f^l ^<?i(^'P^)<-'^ii(^'PK/.). 

where the equalities are obtained by adding terms without residues at the branch points 
to the integrand and using theorem 14.41 One can perform the same trick for the second 
term in Eq. (14-9^ and this proves the lemma. □ 

4.6 Examples. 

Let us present some examples of correlation functions and free energy for low order. 
Correlation functions. 

To leading order, one has the first correlation functions given by: 

W^'\p,q)=B{p,q). (4-44) 



W^°\p,Pl,P2) 



P 



+ P^ 



(4 - 45) 



Res -f^ [B{q,p,)B{q,p2) + B{q,pr)B{q,p2)] 

q^a uj{q) 



Wi°\p,Pl,P2,P3) 



P ► J ► 




+ perm. (1,2,3) 
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+ 




+ perm. (1,2,3) 



1 2 

dEg{p) dEriq) 



Res Res 

q~*a r^a U![q) U[r 



[B{q,pi)B{r, p2)B{r,p3) 



+B{q,pi)B{r,p2)B{r,p3) + B{q,p2)B{r,pi)B{r,p3) 

+B{q,p2)B{r,pi)B{r,p3) + B{q,p3)B{r,p2)B{f,pi) 

+B{q, P3)B{f, P2)B{r; pi)] 
dEg{p) dEr{q) 



+ Res Res 



[B{q,pi)B{r,p2)B{r,p3) 



q^a r^a uj{q) Uj{r) 

+Biq,p,)Bif,p2)B{r,p3) + B{q,p2)B{r,pi)B{f,p3) 
+B{q,p2)B{f,pi)B{r,p3) + B{q,p3)B{r,p2)B{f,pi) 
+B{q,p3)B{r,p2)B{r,pi)] 



(4 - 46) 

First orders for the one point correlation function read: 




Res "^Biqrq) 



(4 - 47) 



q-*a 



uj{q) 



W['\p) = 




Res Res Res ^^^5(r, f)S(., .) 

q,_»a r-*a s->a iiJ{q) UJ{r) i^{s) 

+ Res Res Res ^^^S(r, ^)S(., ^) 

q^a r^a s^a uj[q) Uj[r) Uj[s) 
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+ Res Res Res r)5(., ^) 

g^a r^a s^a LJ[r) UJ[S) 

+B(s, q)B{s, r) + B{s, q)B(s, r)] 



2 P- 



+ 2 P- 



; + P 




where the last expression is obtained using lemma I4.2[ 



Free energy. 

The second order free energy reads 

2F(^) = Res Res Res Res ^M^^^5(r, r)5(., ^) 



(4 - 49) 



p— ►a g— >a r^a s— >a 

+ Res Res Res Res 

p— >a g— >a r— >a s^a 

+ Res Res Res Res 

p— >a g— >a r^a s^a 



c<j(g) uj{r) uj{s) 
^{p)dE,{j>) dEriq) ^p)dEs{r) 



uj{q) uj{r) uj{s) 
<l>{p)dEg{p) dEriq) dEsir) . 



B{r,q)B{s,s) 



uj{q) uj{r) u{s) 
+B{s, q)B{s, r) + B{s, q)B(s, f )] 



B{q,r)B{s,s) 



4.7 Remark: TeichmuUer pant gluings 

Every Riemann surface of genus g with k punctures can be decomposed into 2g + k 
pants whose boundaries are 3g+k closed geodesies (in the metric with constant negative 
curvature). The number of ways (in the combinatorial sense) of gluing 2g + k pants 
by their boundaries is clearly the same as the number of diagrams of Glf \ and each 
diagram corresponds to one pant decomposition. 
Example with k = 1 and g = 2: 



(2) 




5 Variations of the curve 

The goal of this section is to study how the various F^^^ and correlation functions 
change under the variations of moduli of the curve. 

Consider an infinitesimal variation of the curve S ^ S + 6S. It induces a variation 
of the function y{x) at fixed x: 

6fiy\^dx = -fl. (5-1) 
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If we use a local coordinate z, we may prefer to work at fixed z instead of fixed we 
have a Poisson structure: 

b^y\zdx-b^x\^dy = -Vl (5-2) 

The possible fi's can be classified as first type (holomorphic), second type (residue- 
less, and vanishing ^-cycles) and third type (only simple poles and vanishing ^-cycles), 
see [121 for this classification. 



5.1 Rauch variational formula 

Equation 15-21 implies that the variation of position of a branch point is given by: 

b^x[ai) = (5-3) 
dy{ai) 

We assume here that ^ has no pole at branchpoints. Then, Rauch variational formula 
[68| W7\ implies that the change of the Bergmann kernel is 

A M \^ ^(^^) R mr,p)B{r,q) 

^.^(p,g)L(^,,(^) = Lrf^gf.^ dx{r) 

= y Res ^MlMllZMli) . (5-4) 
^ r^ai dx[r)dy[r) 

In particular after integrating over a ^-cycle we have: 

VL{r)B_{r^p)du{r) 



dx{r)dy{r) 
(5-5) 

and integrating again over a B-cycle: 

VL{r)du{r)du^{r) 



^ r^ai dx{r)dy{r) 

(5-6) 

Let us compute the variations of the K-modified Bergmann kernel: 

5nB{p, q) = 5nBip, q) L(p),^(g) + 2m 5^du\p)\^^^^ Kdu{q) 

+2i-ndu\p) K Sndu{q)\^f^g^ 
j^gg ^{r)B{r,p)B{r,q) 
r^a dx{r)dy{r) 

fl{r)B_{r,p)du^{r)Kdu{q) 

+2m Res , . ^ . ^ 

r^a dx[r)dy[r) 

+2' R ^{''^)B_ij'A)du^{p)i^du{r) 

r^a dx{r)dy{r) 
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j^^g ^{r)B{r,p)B{r,q) 

r^a dx{r)dy{r) 

^ p^^^ n{r)B{r,p){Bir,q)-B{r,q)) 

r-^a dx{r)dy{r) 
_^ j^gg ^ir){B{r,p) - B{r,p))B{r,q) 
r^a dx{r)dy{r) 
VL{r-)B{r,p)B{r,q) 

r^a dx{r)dy{r) 

+4 ^ R ^{i^)d'^'^{p)i^du{r)du^{r)K,du{q) 



Res 



dx{r)dy(r) 

^?<!«^-2™*'(p).fer«A(,) (5-7) 

dx{r)dy[r) 



I.e. 



5q + tr (k 6nT K 



B{p,q) 



Res 



x{p),x{q) 



-2 Res 

r—*a 



Q{r)B{r,p)B{r, q) 
dx{r)dy{r) 
n{r)dEr{p)B{r, q) 



(5-8) 



uj{r) 



We thus define the covariant variation: 




(5-9) 



It is more convenient to rewrite Eq. fl5-8p as follows: 

r^a s^r {y{r) — y{r)){x{s) — x{r)) 
= 2 Res Res ^ir-)dErip)B{s, q) 

r^a s^r (?/(r ) — y (r) ) (s(s) — X (r) ) 
fl{r)dEr{p)B(f,q) 



2 Res 

r— >a 



uj{r) 



Res ^^^^ [Q{r)B{r, q) + (](r)S(r, g)] . 



r^a ljj[r] 



(5-10) 



because now we recognize the propagator and vertex of def l4.1[ Similarly, by integrating 
once with respect to g, near a branch point Oj we get: 



dEg{p)\, . 



{p)M'i) 



_2 Res ^^;^n{r)dEg{r) 



Res 



r-*a UJ[r 

dErip) 



r^a iO[r] 



[n{r)dEg{r) + n{r)dEq{i 
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(5-11) 



Those two relations can be depicted: 



D P- 



SI 



p — y 



Q. 



p > — • q 



and 



D P >- 



n 



> ^ > 



From this last variation, one can compute the covariant variations of the correlation 
functions and free energies through the following lemma: 



Lemma 5.1 For any symmetric bilinear form f{q,p) = f{p,(l): 



Dn { Res ^^f{q,q) 
g^a uj[q) 



x{p) 



Res Res — 7^^{r) — f{q,q) 



ij uj{r) 



uj{q) 



(5 - 12) 



a 



D P > <^ 



> * > 



P^ 



Graphically, this means that taking the variation of a diagram just consists in 
adding a leg Q in all possible edges of the graph. In particular if Q can be written as: 



nip)= / B{p,q)Aiq) 
I an 



(5-13) 



where the path dQ does not intersect small circles around branch-pointqj, then we 
have: 



^This excludes the case where il corresponds to the variation of an hard edge, cf [5S1 [TBI [T^ . 
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Theorem 5.1 Variations of correlation functions and free energies: For g + k > 1 we 
have 



^(p«) Jan 



(5-14) 



and, for g > 1, 




(5-15) 

This theorem is proved in appendix [Bl It foUows directely from lemma 14.11 and 
lemma 15.11 

5.2 Loop insertion operator 

In particular for any point q lying away from the branch-points, if we choose 



n{p) = B{p,q) 



(5-16) 



we call DB{.,q) the loop insertion operator, by analogy with matrix models [71 139]. 
It acts on the correlation functions and free energies as follows: 



Theorem 5.2 



and 



Db Wt'\p^, ...,Pk) = W^ii(pi, ...,Pk,q) 
DbF^3) = -wi^\q) 

Db = y{q)dx{q) + — ^— \ k (p ydx] (p (p W^qK (p yd. 

4^7r \ Jq J JbJb ' Jb 



(5-17) 
(5-18) 

(5-19) 



Thus, the loop insertion operator, adds one leg to correlation functions. 



5.3 Variations with respect to the moduH 

Let us consider canonical variations of the curve corresponding to each moduli of the 
curve defined in section 13. 4[ We use Eq. (13-651) 



ydx = ^ ta,kBa,k + ^ tadSa,o + S^TT ^ eidui{p) 

a,k a i 



(5-20) 



to identify the f2's corresponding to varying only one modulus. 
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Variation of filling fractions 

Consider the variation of the curve 

Q{p) = -2indui{p) = - f B{p,q). (5-21) 

i.e. dQ = Bi and A = —1. It is such that: 

6nej = 6i, , 6ntc. = , 6nV^ = 0. (5-22) 
Therefore it is equivalent to varying only the filling fraction = ^ ydx: 

D-2i-Kdu^ = (5-23) 

Theorem 15.11 gives 

A Wi'\p,, . . . ,p,) = - j( Wl^},{pr, ...,pk, q), (5-24) 

and 

<9e,; 



and 

de 



(f Wi^\q), (5-25) 
= - y ydx + ^(^K ^ ydx^ S-2indu, {r) k j ydx. (5-26) 



Variation of temperatures 

Let a and a' be two distinct poles of ydx. Consider the variation of the curve 

pa' 

n{p) = -dSe^^a'ip) = / B{p,q) , i.e. dn = [a, a'] , A = 1 (5-27) 
It is such that: 

5nej = , 5ntf3 = 5a,i3 - 5a',i3 , 5nV/3 = 0. (5-28) 
Therefore it is equivalent to varying only the temperatures ta and ta''. 

Notice that it is impossible to vary only one temperature ta since we have tp = 0. 
Theorem 15.11 gives 



(g) 

k+l 



{pi,. . . ,Pk,q), (5-30) 



and 

^ d 

dtry dt 
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Variation of the moduli of the poles 

Let a be a pole of ydx. Consider the variation of the curve 

= = Res B{p, q)zl{q), (5-33) 

a 

i.e. dVL is a small circle around a and A = tt- -2^. It is such that: 

= , dntfs = , 6ntp^k' = 5a,ph,k' (5-34) 
Therefore it is equivalent to varying only the coefficient ta,k'- 

C'ta,k 

Theorem 15.11 gives 

^ Wt\p,, ...,Pk)= Res 4(g)l^ii(l?i, ...,Pk,q), (5-36) 



^ F^<^^ = - Res z'^{q)wi'\q) (5-37) 



dtr 



and 



_d 

di 



—F^°^ = Res ydxz^ + ( k / ydx] 6b .{T)K(f ydx. (5-38) 
a,k " 427r V / Jb 



5.4 Homogeneity 

Theorem 5.3 For g > 1, we have the homogeneity property: 

(2 - 2^)F(^) = Y.'.,kj^F^'^ + E^"f^^'^ + E^^TT^^'^ (5-39) 

z. e. F^^) zs homogeneous of degree 2 — 2g. 
The proof is given in appendix [Bl 

5.5 Variations of F*^^) with respect to the moduU. 

In this section we compute the ffist and second deriatives of F^^^ with respect to the 
moduli of the curve. This paragraph is only for bookkeeping since those expressions 
have been known for some time [63l [12], [13]. Here we set k = 0. 
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First derivatives of F^^\ 



dta,l3 



— = Kes z^yax 



ydx 



Second derivatives of F^^\ 

52 ^(0) 



{5a,p - 1) Res Res z^{pfB{p,q)zp{q)^ 



piO) 



dta,kdt 



Res z^dS^ 



52 ^(0) 



2i'7T Res z^dui = — B^^k 



dta,kdei a Jig 

52 ^(0) 



52 ^(0) 

deidej 



= 2i'n{ui{(3) - Ui{a)) 
(-2i7rrij /or 



-92 p(0) 

= ln(rfC(«)rfC/3(/?)i?(«,/3)') 
In 



where Co = — is a local coordinate around the pole a. 

Remark 5.1 The definition of given in Eq. (j4-13p is notiiing but tiie homogeneity 
property since it is written in terms of the first derivatives. One can also write a formula 
focusing more on the second order derivatives of F^''-' : 

F^°^ = -J V Res Res q)Vp((i) + V t/3 Res Vo^dSp^o 

a,f3 a,f3 



a,P 

(5-51) 



^ tatp In (7„,/3) - ^e* j( ydx 
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where 

ln7a,a = - / {dSa,o' -\ ^) + In (zq, (o)) (5-52) 

and 

( E{a,P)E{o,o') \ 

ln7a,/3 = In , . . 5-53 

\E{a,o')E{fi,o) J 



One can notice that, in these terms, 

52 ^(0) 

a,fi 



ln(7«,a7/3,/3)- (5-54) 



6 Variations with respect to k and modular trans- 
formations 

6.1 Variations with respect to k 

We have introduced the matrix n in order to easily compute modular transformations 
of our functions. Somehow variations of k, play the role of infinitesimal modular trans- 
formations. Therefore it is important to compute d/dn, and we will use this result in 
section 16. 2[ 

First, notice that W^^\pi, . . . ,pk) is a polynomial in k of degree 3g + 2k — 3, and 
F^^^ is a polynomial in k of degree 3g — 3 for g > I (number of propagators in a graph 

Theorem 6.1 



+ 1 E.ELc;.£.^5.<|Vl(PL,r) 




W^K|Vl(PK/L,^) 






(6-1) 


and in particular for g > 2: 










(r) / W['~''\s) 



(6-2) 

and 

^F« = -. (6-3) 

This theorem is proved in Appendix [B] 

Notice that these equations are to be compared with the Kodaira-Spencer theory 
[Il[2l[IIlll2]. 
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6.2 Modular transformations 

Consider a modular change of cycles: 

where 5^/^ = 5^^,, 5a'b = -^ab'^ ^B'B = 8b'A = -^ba' ^^'^ the matrices 

^AA', ^AB', ^BA', ^BB' have integer coefficients and satisfy 5_aj\,'5qq, — Sba'^ab' — 

Under this transformation of the cycle homology basis, the Abel map and the matrix 
of period change like: 

du' = Jdu , du = J-^du' (6-5) 
with J = {51^, + t'SI^b') = i^BB' - tSab')"^ and 

r' = {Sbb' - rS^>)~^{-SBA' + t^aa') , ^ = {^aa' + ^'^ab')~^ i^BA' + ^'^bb')- 

(6-6) 

Let us define the following symmetric matrix: 

« = = (Sbb'S;^^, - r)-' = 5,^, J (6-7) 

it is such that the Bergmann kernel changes like: 

E = B + 2m du^Kdu. (6-8) 

The K-modified Bergmann kernel changes like: 

-B'(p, q) = B{p, q) + 2m [du'\p)Kdu [q) + du\p)K,du{q)] 

= B{p,q) + 2mdu\p){K + j'nj)du{q). (6-9) 

In other words, the effect of a modular change of cycles is equivalent to a change 
K — > K -I- J^kJ in the definition of the kernel i?(p, q). 

Thus, the modular variations of the free energies satisfy the following theorem: 

Theorem 6.2 For g >2 the modular transformation of the free energies consists 
in changing k to k + J^kJ in the definitions of the modified Bergmann kernel and 

Ahelian differential. 

The first correction F^^^ changes like: 

jrW ^ _ 1 In _ (6-10) 

An equivalent way of saying the same thing, is that: if we change the basis of cycles 
and change the matrix n ^ J*~^{k — 'k)J~^ , then the F^^^ 's are unchanged for g > 1. 
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proof: 

The result for g > 2 comes directly from the variation of the Bergmann kernel. 
F^^^ depends on the cycles only through the Bergmann tau function. Since, 

d\n (tb'x) _ d\n (tbx) _ S^jp,?) - B{p,p) 

dx{ai) dx{ai) p^a, dx{p) 

du\p)Kdu{p) 
= 2in Res -— 

p^ai dx[p) 

du\p)Kdu{p) 
= —ziTT Res -— 

p^a, dx[p) 

— t\i 

dx{ai) 

= -It [Obb' - toab' 



{6-11] 



dx{ai) 
dx{ai) 



and this characterizes the Bergmann tau function totally (up to a general multiplicative 
factor), one obtains the second result of the theorem. □ 

Remark 6.1 The transformation of leading order F^'^^ is more complicated and its compu- 
tation is more involved as the final result depends explicitely on the position of the poles a 
in the fundamental domain. Let us just mention that it depends on all the parameters of the 
modular transformation explicitely. 

Theorem 6.3 // one chooses k, = ^ f— , then F^3\k) is modular invariant. 
proof: 

It is well known that for that value of k, the modified Bergmann kernel is the 
Schiffer kernel and is modular invariant. Indeed, it is easy to check that if k = 
one has k + J^kJ = ^rt^- 

Since the only modular dependence of F^^^ for (7 > 1 is in the Bergmann kernel, 
this proves the modular invariance of the F^^^^s. 

□ 



7 Symplectic invariance 

The following theorem is mostly the reason why we call F^^^^s invariants of the curve. 
This theorem seems to be rather important and it has beautiful applications as we will 
see in section 110.4.11 
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Theorem 7.1 The following transformations of £ leave the F^^^ 's unchanged: 

cx+a ^ ad— be 

• y y + R[x) where R is any rational function. 

• y —>■ y and x —x. 

• y ^ X and x ^ y. 

Notice that these are transformations which conserve the symplectic form 

\dxAdy\ (7-1) 
In particular we have the PSL2{C) invariance: 



x\ fa b\ I X 
y ) \c dj \y 



(ad - hcf = 1 (7-2) 



That symplectic invariance seems to be a very powerful tool to see if different matrix 
models have the same topological expansion. For example, in section 110.4.11 we show 
how symplectic invariance can be used to provide a new and very easy proof of some 
properties of the Kontsevitch integral. 

proof: 

Invariance under the first 2 transformations is obvious from the definitions, because 
the only x and y dependance of the wjf^^s is in uj{q) = {y{q) — y{q))dx{q) which is 
clearly unchanged under the first 2 transformations (notice that the transformation 



X 



conserves the branchpoints). In fact the first 2 transformations leave Wf 



ia) 

k 



unchanged. 

In the 3rd transformation the only thing which changes is the sign of u, and it is 
easy to see that wjf^ is multiplied by (— l)29-2+fc ^ (—1)'', and F^^^ is multiplied by 
(-1)29-2 = 1. 

The 4th transformation is the difficult one. The proof consists in building some 
"mixed correlation functions", and seeing that their definition by inverting the roles 
of X and y lead to the same objetcs. Since it is long and involves new results in the 
framework of matrix model, it is written in a separate paper [13]. 

The case of F^^^ and F*^^^ are done separately in appendix O CH 

8 Singular limits 

Consider a 1-parameter family of algebraic curves: 

£ix,y,t) (8-1) 
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such that the curve at t = has a singular branchpoint a with ap/q rational singularity, 
i.e. in some local coordinate z near a we have: 

t = 

x{z) ~ x{a) + iz- a)i (8-2) 
y{z) - y{a) + {z - a)P 

At t ^ 0, the singularity is smoothed, and we have (the local parameter is now 

xiz,t)r^x{a)+t'^- QiO+o{t'^'^) 

where Q is a polynomial of degree q and P is a polynomial of degree p, and where is 
some exponent which depends on the choice of the parameter t. 
The curve 

^sing(^, r?) = I ^l^j Z p\l] = Resultant(g -^P-v) (8-4) 

is called the singular spectral curve. 

One observes that F^^\S) is singular in the small t limit, and it behaves like 

F(^) {£ (t) ) ~ r« F^^^^g + o(r« ) , for ^ > 2 (8-5) 

F^'\£(t))^-^{p-l)(q-l)u\n(t) + 0(l) , for^ = l (8-6) 

-^s'ing is called the double scaling limit of F^^\ The exponent 'fg and F^^^g are given by 
the following theorem: 

Theorem 8.1 Singular limits: 

) = i^^^^(^^sing) , for^>2 

(8-7) 



and 

^g^i2-2g){p + q)iy. (8-8) 

In other words, our construction of F^^^ commutes with the singular limit. 

proof: 

It is easy to see that the most singular term in the limit of the Bergmann kernel 
in that regime, behaves like and thus dEz'{z) as well. The denominator {{y{z) — 
y{z))dx{z) in the recursion behaves hke i''^*'+^)(P(C) — P{Ci)Q'{C)dC,, and by recursion 
on k and g, we easily see that: 
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if all Zj's are close to a, and is subdominant if some Zj's are not in the vicinity of a. 
The leading contribution to wjf^ is thus obtained by taking z' and l' in the vicinity 
of a only in eq j4-9l i.e. uliI^\(i, . . . ,(k) obey the same recursion formula as eq ]4-9l with 
the curve £^sing- The same holds for the free energy. 
□ 



9 Integrability 

Here, we prove that Z is a tau-function, because it satisfies some Hirota equation. 
9.1 Baker- Akhiezer function 

Given two points ^ and 77 in the fundamental domain, we define the following kernel as 
a formal series in 

-Nf^ydx 

Kn{^,v) = , 



^^p -E E ^N^-''-' I [ ... [ w^'\p,,...,p,) 

\ g=0 l=l,2-2g-l<0 

(9-1) 



where the integration path lies in the fundamental domain. 
This kernel has the following properties: 

• Notice that (a;(^) — x(?7))Kjv(^, 77) — > 1 when t] ^. 

• We have: 

where Wi = ^ N^-^9wi'^\ 

• we have: 

K^{^,r^) = K_^{r^,0 (9-3) 

• One may think that Kn is singular at branchpoints because Ini^'Ar has poles. 



however, using the singular limit theorem 18. ![ we see that the leading behavior 
of all wj'^^^s is given by the Wi^^^s of the Airy curve y = a/x described in section 
110.51 Therefore, near a branchpoint a, when ^,7] ^ a, to leading order is the 
Tracy- Widom kernel [72] : 

J^N[t.,r]) J 



AA 



i = N^/%x{^) - x{a)) , fi = N^/^{x{ri)-x{a)) (9-4) 



In other words, is not singular near branchpoints. 

• The only singularities of K^^C,,!]) are essential singularities at all the poles of 
ydx, with a singular part equal to exp (— J^ydx), as well as a simple pole at 



Then, given a pole a of ydx, we define for ^ in the vicinity of a: 

i^a,N{0 = , (^a(O) 

E{^,a) ^J dx{i)dCa{a) 

-E E In'-''-' / / •••/ wi'\p,,...,p,) 

lini I ./^He-^^'^^"^ (^«(r]))"^*" 1 (9-5) 



rj—ra 



dCaiv) 



where (a is the local parameter near a, Co = j-, and 4>a,N{0 = V^a,-Af(0- 
They have the following properties: 

• '^a,N was defined only in the vicinity of a, but it can be easily analytically con- 
tinued to the whole curve, by choosing an arbitrary point o in the vicinity of a 
and writing: 

/ {ydx - dV^ + t^—) + V;(0 - t« In (z„(0) 

r dZa 

ydx + / {ydx - dV^ + t« — -) + K«(o) - In {za{o)) (9-6) 



• Using the singular limit theorem 18. II near branchpoints, we see that the leading 
behavior of all VT/^^'s is given by the VT/^'^'s of the Airy curve y = \fx described 
in section 110.51 Therefore, near a branchpoint a, when ^ a, we have: 

i;^,N{0 C At{i) , i = N'/'{x{0-x{a)) (9-7) 

where C is some normalization constant {C = ipa,N{O')/A'i{0))- In other words, 
ipa,N is not singular near branchpoints. 

• The only singularities of 4'a,N are essential singularities at all the poles of ydx, 
with a singular part equal to exp (— ydx). 

This is why we call those formal functions Baker- Akhiezer functions (cf [9]). 
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Remark 9.1 In fact, those functions are exactly Baker- Akhiezer functions only when the 
curve has genus ^ = 0. In the general case, the Baker-Akhiezer functions must also have 
the property that they take the same value after going around a non-trivial cycle. It is 
not difficult to multiply il)a,N by the appropriate 6 function in order to fulfill that property. 
However, if we do that, we destroy the l/N"^ expansion. 

This is why the ipa,N defined above can be called a "formal Baker-Akhiezer function". 
This definition is sufficient to find a formal Hirota equation, valid only order by order in 

9.2 Sato relation 

Given two points ^ and t] on E, and a complex number r, we define the curve: 

£ + r{i,-r,] = \J^x{p)Mp) +r'^^^^ , p G s| (9-8) 

The differential ydx + rdS^^^j has the same v4-contour integrals as ydx, the same poles 
with the same singular part, plus two additional simple poles, one located at p = ^ 
with residue r, and one at p = r] with residue — r. 
We have Sato's relation: 



Theorem 9.1 



kn{^,v)= ' ^ ;;r , mo= ' ^ " (9-9) 

Indeed, the definition of Kj^ is the formal Taylor expansion in powers of r = of 
the RHS. 

9.3 Hirota equation 

Consider two algebraic curve S{x,y) and S{x,y) with the same conformal structure 
(i.e. the same compact Riemann surface S), then we have an Hirota bilinear relation: 

Theorem 9.2 We have the bilinear relation: 



and also: 



Res dx{r]) K^i^, r]) Kf,{r,, Q = K^{i, Q (9-10) 



Res dx{i) ^«,jv(0 ^„ = if Nt^ >Nt^ + l (9-11) 



which takes exactly the form of the Hirota equation /3, \6U^ . 

It is important to notice that this Hirota equation makes sense only order by order 
in its l/N"^ expansion. In the way we have obtained it, it is meaningless for finite 
(appart maybe from the genus zero case ^ =0, under the condition that the 1/A^^ 
series is convergent). Therefore, we have a "formal Hirota equation". 
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10 Application: topological expansion of matrix 
models 



In this section, we show how the objects defined in section J4.2.4l (i.e. k = 0) correspond 
to the terms of the topological expansion of the free energy and correlation functions 
of various matrix models when one considers appropriate curves S{x,y). Notice that 
in all this section we consider k = 0. 



10.1 Formal 1-matrix model 

The formal 1-matrix model (cf [H]), is known to be the generating function which enu- 
merates maps of given topology since the work of [Q] , then [H [25l [52] . Its topological 
expansion was computed in several steps. The authors of [7] introduced a recursive 
algorithm to compute the F'-^^'s, only in the 1-cut case (i.e. ^ =0), and then the 
method was extended to other cases [3lll]. The computation of the subleading term F^^^ 
was done in general by Chekhov [15]. The computation to all orders of the correlation 
functions was found in [33j, and the free energies in [T7j. 

10.1.1 Definition 

Definition 10.1 Formal 1-matrix model. 

Consider a "semi- classical" potential V{x), i.e. such that V'{x) is a rational func- 
tion. Let D{x) be its denominator, i.e. D{x)V'{x) is a polynomial of degree d, and let 
Xi, . . . , Xd he its zeroes: 

d 

D{x)V\x)=\{{x-X,) (10-1) 

We write: 

5V,{x) = V{x) - V{Xi) - U^"{Xi){x - X,f (10-2) 
Choose an integer n, and a d— partition of n, n = {ni, . . . , n^}, such that 

d 

J]n, = n. (10-3) 

i=i 

The following gaussian integral (where each matrix Mi is of size rii) is a polynomial 
in T of the form: 



JJ det(Mi ® In, - In, ® Mjf (J2 SV,{M,)) 



^ i=i 



i>j 



47 



dl/2 

= 5]Afc,,T'= (10-4) 

k=l/2 

We define the formal matrix integral as the formal power series in T : 

oo 2k 

Z^MM = Y,T\Y,Ak,j). (10-5) 

fc=0 j=0 

One can also define its formal logarithm, i.e. the free energy 

oo 

-^iMM = — In ^iMM = ^ ^'^-Bfe (10-6) 

fc=0 

It is a standard computation discovered by 't Hooft (f71\\41f ), that, for fixed ei = 
for every k, n~'^Bk is a polynomial in 

^T^^^ /n\ 2-23 

5fc(ni,...,nrf) = 2^ Bk,g{ti,. . .,td) y-j (10-7) 

9=0 

Thus we define the following formal power series in T: 

oo 
k=0 

Remark: Here, the question of convergence of those series is not relevant. It is 
well known that each -Fimm ^ convergent series (because it is written in terms of 
algebraic functions of T), but -Fimm is not. 

10.1.2 Loop equations and classical spectral curve 

It is easy to see (this property holds for each power of T because it holds for gaussian 
integrals) that the formal matrix integral satisfies, order by order in powers of T, the 
loop equations (i.e. Virasoro constraints), which can be written (see [261 [23]): 

y{xf + — X, x) = —j^ ( Tr ^ ^ / ' ) 10-9 

4 n \ X — M / 

where y{x) = - ^(Tr^^) and uj2{x,x') = ( Tr Tr ^7^)^. where the 

expectation value < . > is defined in a formal way similar to F. 

If one identifies the coefficients of n° in each side, one gets an algebraic equation 
(here hyperelliptical), which is called the "classical spectral curve": 

£imm{x, y) = D{xf (y2 _ ^V'\x) + P{x)) (10-10) 
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where D{x)P{x) is a polynomial of degree at most deg{D{x)V'{x)) — 1, and completely 
determined by the condition that the polynomial P{x) is a formal power series in powers 
of T such that at T = 0: 

P{x,T = 0) = J2^,^^ (10-11) 
It is such that there exist some contours Ai, i = 1, . . . ,d such that: 

^ / ydx = e, (10-12) 

Notice that the genus ^ of the curve £^imm is the number of non vanishing e^'s minus 
one. 

Most often in the litterature, V is chosen polynomial such that V'{0) = 0, and only 
the 1-cut case is considered, with only one non-vanishing filling fraction at X = 0. The 
resulting curve has genus ^ =0. This is the case which is relevant for enumerating 
polygonal surfaces. 

It was proved in [331 [U] ) the case of polynomial potentials only (but it is clear 
that it can be extended to the semiclassical case), that one has: 

Theorem 10.1 

I (10-13) 

This proves that F^^'' (which we recall is a formal power series in T) genericaly has 
a finite radius of convergence T < T^. 
We also have: 



Tr . . . Tr \ = \^ N' 



x{pi)-M"' x{pk) - M / ^ dx{pi) . . . dx{pk) 

+h,iNi^V\xip^))-yipi)) 

5k,2 



(10-14) 

10.2 2-matrix model 



(x(pi) - X{p2)y 



The formal 2-matrix model is known to be the generating function which counts bicol- 
ored maps (let us say the 2 colors are -|- or -, thus it is an Ising model on a random 
map), it was introduced by Kazakov [51]. The loop equations were first written in [69] . 
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F^'^^ was computed in [631 [IS [IS [38] . F*^^-* was first found in pi] for the ^ =0 case, 
the in [35] for ^ =1, then in [37] for arbitrary ^ . Then higher orders for the corre- 
lation functions were first derived in [ID], and the F^^^^s for g > 2 were first found in 
|19] . During the same time it became clear that matrix models topological expansion 
was closely related to algebraic geometry [HSl [221 [23, [21] ■ 

10.2.1 Definition 

Definition 10.2 Similarly, consider V( and V2 two rational functions with respective 
denominators Di{x) and D2{x). The equation: 



has d = deg{V(Di) * deg(V2-D2) solutions. We then write: 

SViA^) = Vii^) - yliX^) - Yi{x - Xi) - ^^^^^{x - Xi)^ (10-16) 

et 



V"(Y) 

SV2M = V^iy) - V^iY^ - X,iy - Y,) - ^f^iv - Y^)' (10-17) 
We then choose an integer n, and a d— partition of n: 

d 

n = ^ni. (10-18) 

i=l 

The following gaussian integral (where each matrix Mi or Mi is of size rii) is a 
polynomial in T of the form: 

t})^ e-f Tr mx,)+v,iYO-x,Y,) f ^M, . . . dM.dM, ...dM, 
/!T' J 

iTr (^ll^(M,-xanj2+l2^(Af,-yanj2-(Af,--xan,)(A/>-yan,)) 



i=l 

i>j i>j 

TT6V,,,{Mi) + 6V2AM^)y 



dl/2 

= J]Afc,;T'= (10-19) 

k=l/2 

Similarly to the 1-matrix case, we can define the formal 2-matrix model as a formal 
power series in powers of T (see l4^ ): 

00 2k 

(10-20) 

k=o j=o 
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One can also define its formal logarithm, i.e. the free energy 

oo 

F2MM = — In Z2MM = ^ T'^Bk (10-21) 



fe=0 

Again, it is a standard computation (l71[ \4^), that, for fixed ti = for every k, 
n~'^Bk is a polynomial in \jv?: 

5max(fc) „ „ 

5fc(ni,...,n,)= Bk,gieu...,ed) [-) (10-22) 

9=0 

Thus we define the following formal power series in T: 

00 

F^MM = Y.^'B,,9{^i.---.^<i) (10-23) 

k=0 

10.2.2 Loop equations and classical spectral curve 

Again, the formal 2-matrix integral satisfies, order by order in powers of T, the loop 
equations (i.e. W^-algebra constraints), which can be written (see [2S1I21]): 

{y - y{x))U{x, y) + ^f/(x, y- x) = E{x, y) (10-24) 



1 Vi{y)-Vi{M2) 



where y{x) = Vl{x) - ^ ^ Tr U{x, y) = x - Vi{y) + I ^ Tr ^ 

and U{x,y-x') = ( Tr ^ ^^Mz^ Tr ^)^, and E{x,y) = {y - Vl{x)){x - 
^2 (Z/)) - - ( Tr v^{y)-vi{M2) \ ^ j.^ ^^^^^ expectation value < . > 



y-M2 



x-M-i y-M2 

is defined in a formal way similar to F. 

If one chooses y = y{x) and identifies the coefficients of in each side, one gets 
an algebraic equation S2MM{x,y) = 0, which is called the "classical spectral curve" 



S2MM{x,y) = D,{x) D2{y) ((K(x) - y){V^{y) - x) - P{x,y) + T) (10-25) 

where -Di(x) D2{y)P{x, y) is a polynomial of degree < deg{DiV() in x and a polynomial 
of degree < deg{D2V2) in y, and with fixed filling fractions: 

i-£,dx = r^ = ., . ^., = r (10-26) 

and such that in the limit T ^ and V/ e/ 0, one has: 

y-V;{x)-J2-^ + 0{T') , x^V^{y)-Y.7^.+0{T'). (10-27) 
X -A-i . y ±1 

Then one has: 
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Theorem 10.2 



2MM 



2UM 



(10-28) 



proof: 

This theorem was proved in [ICTl [T^ . □ 

Again, this proves a posteriori, that -F2MM (which we recall is a formal power series 
in T) genericaly has a finite radius of convergence T < T^. 
We also have: 



Tr 



1 



x{pi) - Ml 



Tr 



x{pk) - Ml 



TT/(9) 

T2-k-2g nLk 



9=0 



w:ripi,---,pfc) 



(ix(pi) . . . dx{pk) 
+5fc,iiV(\/'(x(pi))-y(pi)) 

4,2 



;i0-29) 



(x(pi) - x{p2)y 



Remark that the 1-Matrix model is a special case of the 2-matrix model with V2 a 
quadratic polynomial. 



10.3 Double scaling limits of matrix models, minimal CFT 

It is well known that double scaling limits of the 1-matrix model, or 2-matrix models 
are in relationship with (p, q) minimal models of conformal field theory [59l [26| [20] . 

We have seen that as long as the curve is reg ular, all the F^^^'s can be computed. 
This shows that the radius of convergence in T of F^^\T) is reached for singular curves. 
So far, only rational singularities have been studied in detail. 

Thus, consider the case where the potentials Vi and V2 are fine-tuned so that the 
curve S2MM has a. p/q singularity at T = Tc (notice that for the one matrix model one 
necessarily has g = 2): 

r T = T, 

I x{z) ^p^a x{a) + {z - z{a))'^ (10-30) 
[ y{z) r^p^a y{a) + z{a))P 

We can use the notation of section [8] with t = T — T^, and thus at t 7^ 0, the singularity 
is resolved, and we have (the local parameter is now ( = zt^^): 

x{z,t) ~ x(a) +t«'^ QiO + oif'") 

y{z, t) ~ y{a) + t^^ P(C) + oit^ (10-31) 
C = zt-" 

where Q is a polynomial of degree q and P is a polynomial of degree p. 
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The curve 



^sing(e, ^) = { ^ = = Resultant(Q -^,P-v) (10-32) 



is called the singular spectral curve. 

The dependence on T of the 2-Matrix model is such that: 



dSoo.,ooy (10-33) 



where oOa; and oo^ are the 2 common poles of x and y. They are far away from branch- 
points, and in particular from the singularity. This means that in the vicinity of the 
singularity we have: 

Ct^dC + Oit^") (10-34) 



where C is some constant o: 
implies the Poisson relatiorl^ 



order 1. After substitution with the limit eq JlO-311 this 



pPiOQ'iC) - qQiOP'iC) = ^ t'-('^^'^-'> (10-35) 

and therefore 

u = J (10-36) 
p + q -1 

Thereom 18.11 proves that: 

^2Mm(^) ~ (T-T,)(2~2.)(p+,)/(p+.-i)^(.)(^^.^j ,for^7>2 (10-37) 

We also have: 

f'2MM(r) ^-^^ ~ (P - 1)(9 - ^> 'n (r - + 0(1) , fors = 1 (10-39) 

We thus have a way to compute explicitely the double scaling limit of -F2mm- 

10.3.1 (p,q) minimal models 

Let us study in more details the curve £{p,q) (cf PU]). 

As we have seen above, the curve for the (p, q) minimal model is of the form: 

£(p,,){x, 2/) = { ^ I p(^) = Resultant(Q -x,P-y) (10-40) 



^This Poisson relation is well known and can be found in 
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where P and Q are polynomials of respective degrees p and q, satisfying: 

pPQ' - qQP' = ^ (10-41) 
The solution of which can be written |26j : 

P = (QP/'?)+ (10-42) 

and 

(QP/</)_ = ^^1-9 + OiC'), (10-43) 

q 

where we have used the notations f = f+ + f-, with /+ and /_ denoting respectively 
the positive and the negative part of the Laurent series of /. This last equation implies 
g — 2 equations for the coefficients of Q. 

The curve £^{p,g) has genus zero ^ =0, and is such that x and y have only one pole 
a = oo. The Bergmann kernel is: 

B{CuC.) = jP^,. (10-44) 

(Ci - (2)^ 

The moduli (of the pole) of that curve are the Qk and P^ such that: 

QiC) = i2Q^^' ' nC) = X;i^.C' (10-45) 

by a translation on (, we can assume that Qq^i = 0, and by a rescaling of ( we can 
assume that Qq-2 = —qQq, and the Poisson Eq. (I10-4ip implies that Pp-i = and 
Pp-2 = -pPp, thus: 

Qq-i = P,-i = o , %^ = -g , ^ = -P- (10-46) 

We find: 

F(°)(£:(,,,)) = 0. (10-47) 

10.3.2 Other times 

More generaly |26|, we can deform the {p,q) minimal model with p + q — 2 times 
ti, . . . ,tp+q_2. For this purpose, one considers Q{C) a degree q monic polynomial, 

q-2 

Q(C) = C + 5^ (10-48) 
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whose coefficients U2, ■ ■ ■ ,Ug are determined as functions of q—1 parameters, ti, . . . , tg-i, 
by the following requirement: 

(gp/'?)_ = y .g-^/'? + + 0{C')- (10-49) 
Then we define the degree p monic polynomial -P(C) by: 

^ = + E ^p-^<' = ^+ ' - E ^^.+.-1^+ " (10-50) 

j=0 i=l ^ 

which depends on p — 1 other times tg, . . . , tq+p_2. 
The corresponding classical spectral curve is: 

£^(p_q)(x, y) = Resultant(x — Q,y — P). (10-51) 

and depends on times ti, . . . , One can check that if ^2 = ^3 = • • • = tp+q-i = 0, 

one recovers the {p, q) minimal model. 

It is well knwon that this curve is the spectral curve of the dispersionless Witham 
hierarchy [6^ . 



10.3.3 Examples of minimal models 

• Airy curve (2, 1) 

The classical spectral curve for the (2, 1) minimal model is 

£i2,i){x,y)=y^ + h-x. (10-52) 

g(C)=C' + ti , nC) = C (10-53) 

It is studied with particular care in section 110.51 since it describes the behaviour of 
a generic curve around the branch points, and thus coincides with Tracy- Widom law 



• Pure gravity (3,2) 

g(C)=C'-2t; , P(C)=C'-3< , h = 'iv^ (10-54) 
The classical spectral curve is: 

^(3 2) {x,y)=x^ - ^v'^x -y^ + 2v^ (10-55) 
and is studied in details in section 110.61 . 
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• Ising model (4, 3) 

Q(C) = C' - 3< - 3«; , P(C) =C^-4<2_4^^ + 2t;2-^t5((2_2t;) (10-56) 
with: 

= 4v^ + Qw^ , t2 = Qvw. (10-57) 
The classical spectral curve is 

%3)(2;,1/) = - y3 - 4t; V + 3 A + 2^;' 

+12tt;t>(— a;y + v x) + Qw {—x + 2vy — Av ) + 8w x — 3w . 

+5t^{—x'^y — v^x^ + 2v'^y + 2v^ — 2wyx + 2v'^wx + Sw'^y — 17v'^w'^) 
25 

H tliv'^y + 2v^ - Avwx - Uvw"^) 

3 

125 

+— t^(-a;2 + 2v^ - 6wx - 9w^) (10-58) 

The variations with respect to the moduli t^, t2 and ti correspond respectively to 
the variations of the form ydx 

= -rfAg = d{C^ - 5vC^ + lOv^C), (10-59) 

= -rfAi = diC + ^ i2v'C + 2vw - wC')), (10-60) 

12 f — 

and 

= -dA^ = d{C + I {v'e - ^vwC - 2w^)) (10-61) 

in the notations of section O 

• Unitary models (g + 1, q) 

S(^q+i,q){x, y) = Tg+i(x) - Tq{y) (10-62) 

g(C)=T,(C) , P(C) = T,+i(C) (10-63) 
where Tp is the p**^ Tchebychev's polynomial. 

10.4 Matrix model with external field 

Define the formal matrix model with external field [781 : 

Zuext 



J rfMe-^T^(^W-^'^^) (10-64) 
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where V'{x) is a rational function with denominator -D(x), and A is a fixed N x N 
matrix. Consider its topological expansion (in the sense of a formal integral as in the 
1-matrix case above [41]): 



oo 



FMext = - In Zue.t = ^Mexf (10-65) 

9=0 

Let us assume that A has s distinct eigenvalues Ai, . . . , As with respective multi- 
plicities mi, . . . , rris such that rrii = N. The minimal polynomial of A is: 

S{y) = ll{y-X.). (10-66) 

i 

Define the classical curve obtained by "removing" the connected term in the 

loop equations: 

SueA^, y) = {{V'{x) - y)S{y) - P(x, y))D{x) (10-67) 



where 



P(,,„ l/T.IM^i;M^M^) (10-68) 

^^^A^\ x-M y-A' ^ ^ 



so that P{x,y)D{x) is a polynomial in both x and y. 
Then one has: 

Theorem 10.3 



41 = Z^^H^Mext) ^ ^ 

;i0-69) 
proof: 

This theorem is proved in appendix [Dl and the proof is very similar to that of the 
2-matrix case above, see 1191. □ 



10.4.1 Application to Kontsevitch integral 

The Kontsevitch integral is known to be the generating function which computes in- 
tersection numbers of moduli space of Riemann surfaces (see [58| [21]). It is defined 
as: 



^ Kontsevitch 



[ dM e-^T^ (^-MiK^+t.)) = 1 Tr 1 (10-70) 

J N K ^ ' 

where A has eigenvalues Ai, . . . , Aat. Thus it is .^Mext with V{x) = x'^/3 and A = A^+ti. 
Its classical curve is: 

^Kontsevitch (a;,?/) = (x^ - y)S{y) - xSi{y) - S2{y) (10-71) 
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where Si{y) and 5*2 (y) are polynomials in y of degree at most s — 1. 

If we assume that the curve has genus zero (which is the case if we want the 
-^Kontsevitch ^e the generating functions for intersection numbers), then we can find 
explicitely a rational parametrization: 

i^Kontsevitch(a:, = | ^ ^2 rj^ ^ "^"^ . (10-72) 



Using the symplectic invariance of theorem 17. 11 we may exchange the roles of x and y. 
There is a unique branch point in y, solution of y'{z) = 0, located at 2; = 0. Since the 
formulae for F^^^ consist in taking residues of rational functions at the branch point, 
we may consider the Taylor expansion of x{z) near z = 0, i.e. 

^:Kontsevitch(x, v) = { Z ^^'^i^'^^" (lO'^^) 

where we have defined the Kontsevitch times: 

4 = lTrA-'= (10-74) 



so that: 



-^Kontsevitch -^'■^H^Kontsovitch) 



(10-75) 



Again, using symplectic invariance of theorem 17. 11 we may add to x{z) any rational 
function of y{z), i.e. we immediately get a 1-line proof of the following well-known 
theorem: 

Theorem 10.4 -^Kontsevitch dtpeuds only on the odd times t2k+i, with k < 3g — 2: 



-^Kontsevitch ~ Kontsevitch I'' 1 ' ''35 '^5, ■ ■ ■ , ''6g-3) 



[10-76) 



And, if we assume that = for k > p + 2, the curve is: 

x(z) = Z — \ 'Y7k=Q'^k+2Z^ 



(10-77) 



y{z) = z' + ti 

which is identical to the curve of the (p, 2) model of section 110.3. H and which is well 
known to satisfy KdV hierarchy. Thus, again we have a 1-line proof of the well known 
result: 

Theorem 10.5 .^Kontsevitch is tt KdV hierarchy tau function. 
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• Examples: the first few correlation functions 

For Kontsevitch's curve we have: 

B[z,z') = j^^^ , dE,{z') = , uj{z) = 2z^dz{2 - ^t^.+s^'') 

\ z z \ z z 

(10-78) 

and the only branchpoint is located at ^ = 0. 

From the def J4.2l we easily get the first correlation functions: 



(10 - 79) 



Wl'\z,,z^,z^) 



(10 - 80) 



dzi dz2 dzs 

^2 ^2 
^1^2 ^3 



Wi'\z,,Z2) 



dzi dz2 



{2-t^)\^zl + bzl + ^zlzl'' 



l{2-hYzlzl 

+Qtlz^z^ + (2 - h){Qhz^zl + Qt^zlzl + bt-jzlzl) 



(10-81) 



and 



W[^\z 



dz 



— -r 252t^;z« + 12t^;z''(2 - t^)mtjz' + 21 ts) 

128(2 -t3)V0L ^ 5 V 3A 7 5) 

+z\2 - t3)2(252t^ + 348 ^5^7^^ + 145 1?/ + 308htgz^) 
+2^(2 -ta) (203^5 + 145 z% + 105 z% + 105 z%i) 
+105 (2 - ts)^^ 



(10-82) 



The first and second order free energies are found: 



4oLevitch = -^lnfl-|) (10-83) 



24 V 2 
and 

^(2) ^ J_ 252ti + 435t5t7(2-t3) + 175t9(2-t3)^ 

Kontsevitch ^929 (2 - ^3)^ ' I ^ O J 

which coincide with expressions previously found in the litterature 
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10.5 Example: Airy curve 

The curve y = ^/x is particularly important, because it corresponds to the leading 
behaviour of any generic curve near its branch points. It is also the minimal model 
(1,2) (cf [261EH]), also called Tracy- Widom law [72]. 

Consider the curve 

£{x,y) = y^ -X. (10-85) 
We chose the uniformization p = y: 

^f^l = (10-86) 

y{p) = p 

There is only one pole a = oo, and there is only one branch point located at a = 0, 
the conjugated point is p = —p. The Bergmann kernel is the Bergmann kernel of the 
Riemann sphere: 

It is easy to see that all correlation functions with 2g + k > 3 are of the form: 

W^'\pi, ...,pk)= 4'\pI ■ ■ ■,pI) dpi... dpk (10-88) 

Moreover, the diagrammatic rules are clearly homogenous, so that the function W^^^ (p) 
must be a homogeneous function of p. It is easy to find that: 



<)(p) = ^ (10-89) 



and the total 1— point function is 



Wi{p,N) = -Nydx + ^N^'^mi^\p) = WiiN^^p, 1). (10-90) 



Similarly, the total 2— point function is 



W2{p,q,N) = ^A^~^Wi^^(p,g) = W2{N^ip, N^^q,l) (10-91) 

g=0 



and in general 



Wkipi, . . . ,Pk,N) = J2N'~''~'wli'\Pi^ ■ ■ ■ ^Pk) = WkiN'^p,, . . .,mpk, 1). (10-92) 

9=0 

The solution of the recursion def J4.2l can be found explicitely in terms of the Airy 
function. 
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Consider g{x) = Ai'{x)/Ai{x) where Ai{x) is the Airy function, i.e. g'{x) +g^{x) = 
X = p^. In terms of the variable p we write: 

f{p) = g{p') , f + f = p'. (10-93) 

2p 



It can be expanded for large p: 



32p 



°° 1 q 



A:=0 

where the coefficients in the expansion satisfy 

4 3 ^ 

/fc-i + f^f>^-^ = '^'^'0- (10-95) 

3=0 

The solution of the recursion def l4.2l for the 1— point function is: 
w( ^^ y - f(P)f(-P) , 2, ^ dp , mdp , 15!! 

/IN . UXp) - fiP'))ifi-p) - fi-p')) , I nnn7\ 

l^.fep , 1) - -4 _ ^ ^^^^^ _ Pdpp dp (10-97) 



In particular 



so that: 



^^(P>P>1)= (/J)'!^7(-P))^ ^^' (10-98) 



1) + 1)^ = V = a; rfa;2 (10-99) 
Similarly we find for instance (with obvious cyclic conventions for the indices): 

^"3(^1,^2,^3, 1) 

dxidx2dx3 

{pI - pDipl - pDipl - pI) 

Ell fiP^)fi-P^)ifiP^-l) + f {-P^-l) " /(Pi+l) " /(-p^+l)) 
(/(Pl) - /(-Pl))(/(P2) - /(-P2))(/(P3) - /(-P3)) 

_ dpi dp2 dp3 ^ dpi dp2 dps ^ 
2pIpIpI 2^pIpIpI 

(10 - 100) 

and one can easily find similar expressions for all W^s. 

In fact all correlation functions can be written with a determinantal formula [3ll [32] , 
with the Tracy- Widom kernel [72] : 

OC Ob 
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The fact that the Baker-Akhiezer function is Ai{x) and satisfies the differential equa- 
tion Ai" = xAi can be seen as a consequence of the Hirota equation theorem j9.2[ 



Remark 10.1 To large N leading order the first term can be written in terms of 

Ferrer diagrams (Young diagrams): 

j O |A|=fc-3 j ^' '' 

where A is a Ferrer diagram, nj(A) = / \j = i}, and Mx are the elementary monomial 
symetric polynomials: 

For instance with k = A there is only one diagram (1), and M(i)(zi, Z2, z^, z^) = zi + Z2 + 
Z3 + Zi, and: 

„.(0), s 3 dpidp2dp3dp4 / I 1 1 1\ nn^nA\ 
Wl '{pi,P2,P3,P4) = 7 2 2 2 2 — + — + — + —• 10-104 

4 Pipipipi \pi Pi Pi pi 



Remark 10.2 The free energies are all vanishing for that curve: 

F^9) = 0. (10-105) 

10.6 Example: pure gravity (3,2) 

In this section we study in details the (3, 2) minimal model, also called pure gravity 
It corresponds to the curve 

{x{z) = — 2v 
y{z) = z^-^vz . (10-106) 

We recognize Tchebychev's polynomials T2 and T3, which satisfy the Poisson relation 
Eq. fllO-4ip . Up to a rescaling z = y/vp, the curve reads: 

{x{p) = vip^ — 2) 
y{z)=v^'\p' -?>p) . (10-107) 
h = 

There is only one x-branch point at p = 0, and the conjugated point is p = —p. The 
Bergmann kernel is the Bergmann kernel of the Riemann sphere: 

62 



^(q) = {y{q) - y{q))dx{q) = iv^^^ (g' - 3) dq 



Under a variation of ti we have: 



dy{p)dx{p) 



dti 



,1/2 



dp = — 1>^/^ Res qB{p, q) 



Ai(g) 



so the effect of d/dti is equivalent to: 



d 



(g) 



10.6.1 Some correlation functions 

Using def J4.2t we find: 

W^3°^(Pi,P2,P3) = 
W['\p) = - 



Vl'^ Res qW^ 

q~*oo 



(9) 
k+1 



(10-109) 
(10-110) 

-v'/'q 
(10-111) 

(10-112) 



V dpidp2dp^ 



6 pIpIpI 



(12)2 p4 



dp 



25 33 p6 qG 



15g^ + 15p^ + 6pY + "^pY + 9p^g^ + Gp'^q^ 

,,A2)i . -15/2 ^ 135 + 87p^ + 36/ + 12p6 + 4/ ^ 
W[\p) = -v I 7 ^1035^10 ^'P 



W^i°^(Pl,P2,P3,P4) = 2^2 2 2 I ^ + ^ 5Z ~2 ) dpidp2dp3dpi 

yPlP2P3P4 \ A Pi 



dpdq 



(10-113) 
(10-114) 

(10-115) 
(10-116) 
(10-117) 



(Pi,P2,P3,P4,P5) _ v-'"/' fi + 3y l + eV^ + sV 1^ 

dpidp2dp3dpidp5 QpIpIpIpIpI \ ^ Pi Pfp] i Pi J 



(10-118) 



etc... 



Using def J4.3l and eq JlO-112] we find from eq llO-1131 

53^(0) 1 1 ^2^(0) ^1/2 



dtl 



and using eq JlO-11^ 



Qv 2v^ dtl 73 

1 92^(1) 1 



1 



dF^_ ^ 

dtl ~ (12)2t.2 ~ 48ti ' dtl ~48t2 

63 



(10-119) 



(10-120) 



as well as using eq JlO-1161 

dF('^ _ -7 7 _ 7 d^F(^^ _ 49 

9ti 2835 28 33/2tp ^ dtj 2933/2^^/2' ^ ^ 

We may thus verify that the second derivative of the free energy: 



satisfies the Painleve equation to the first orders: 



+ = W (10-123) 
D 3 



It is well known that this equation is satisfied to all orders [26], and here this can be 
seen as a consequence of the Hirota equation theorem. 19.21 



11 Conclusion 

In this paper, we have constructed an infinite family of invariants of algebraic curves. 
By construction, these invariants coincide with the topological expansion of matrix 
integrals in the special case where the algebraic curve is the large part of the matrix 
integral's spectral curve. But we emphasize again that the construction presented here 
goes beyond matrix models. 

Our invariants are defined only in terms of algebraic geometry, and they have many 
interesting properties, like homogeneity, and integrability (they obey some Hirota equa- 
tion). 

The problem of computing the F^^^'s for matrix models is an old problem which 
was addressed many times, and which found more and more elaborate answers We 
claim that ours is more efficient, because it contains all multicut cases, and various 
types of matrix models at once. Also, even in the simplest cases (1-matrix model, 1 
cut), our expressions are simpler than what existed before [7j. Our F^^^^s are defined 
recursively, like those of [7], but the recursions are much easier to handle, and it is 
much easier to deduce properties to any order g from our construction. 

The efficiency of our method becomes striking when one wants to compare different 
models (Kontsevitch and KdV for instance), or when one wants to take singular limits. 
Another important application of our method, is to prove that our F'^^^'s provide a 
solution to the holomorphic anomaly equations of [11] in topological string theory, 
thus confirming the Dijkgraaf Vafa correspondence. We claim that this can be proved 
easily from our work, and we present it in a coming paper [12]. It would also be 
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interesting to compare our free energies with the D-Modules considered in [5l [6] as 
partition functions of a unified matrix M-theory by checking that they indeed satisfy 
the equations of [6]. 

One of the reasons our method is very efficient also, is that it can be represented 
diagrammatically, without equations, and thus very easy to remember. 

Perspectives and generalizations: 

• The first thing one could think about is to understand what our F^^^^s compute 
in algebraic geometry. There has been some attempts to recognize the first few 
of them F'-^^ as the Dedekind function [29l[30], F^^-* as the Eisenstein series [50] , 
but the answer for higher g is still obscure. Also a combinatoric interpretation is 
missing. 

• Beyond that, it would be important to understand the combinatorics behind our 
diagrammatic construction. Since all diagrams are obtained from trees, it seems 
to be related to Schaeffer's method for counting maps [70], but this issue needs 
to be investigated further. 

• Double scaling limits of matrix models are in relationship with conformal field 
theory (CFT), and in particular minimal models. It would be interesting to 
compare our formulae with those obtained directely from CFT [281 [ZS [ZZ] , and, 
since higher genus CFT is far less known, maybe our method can bring something 
new to CFT. 

• Also, the link between formal matrix models (i.e. those defined as combinatorial 
generating functions, for which it makes sense to consider a topological expansion, 
cf. [H]), and actual convergent matrix integrals needs to be better understood. 
The difficulty lies beyond any order in perturbation, and integrability could play 
a key role in understanding the relationship in more details. 

• It would be interesting to check that the topological expansion for the chain of 
matrices [38], is also given by the same F^9)^s. This is an open question, but 
there are strong evidences that the answer is positive. For instance it is easy to 
see from [38] that and some of the first few correlation functions are indeed 
the same. 

• It would be interesting also to extend our construction to other types of matrix 
models, for instance non-hermitian (real symmetric or quaternionic). A first 
attempt was done in [18] . Another possible extension is towards the O {n) matrix 
model, whose large limit is known in terms of an agebraic curve [HJ |15] . To 
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leading order, the 0{n) matrix model solution looks very similar to that of the 
1-matrix model, except that correlation functions are no longer meromorphic 
functions on the curve. Instead they gain a phase shift after going around a 
non-trivial cycle. This would probably allow to define a twisted version of our 
construction. 
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A Properties of correlation functions 

In this section we prove the theorems stated in section 14. 4[ 
We use very much the following obvious properties: 



{x{p)-x{q)y 

dE^jp) _ dE,{p) 
uj{q) uj{q) 

and 

dE,{p) ^ , ^ 1 . B{p,q) 



(A-2) 



lim — —— dx{q) = — lim dx{q). (A-3) 

q^ai uj[q) 2 q^ai dy[q) 



The third one is nothing but De L'Hopital's rule. 

Theorem 14.11 proof: 

From the definition 14. 2[ we have: 



W^f(p,Pi,P2) = Res^^(5(g,pi)i?(g,p2)+i?(g,P2)i?(g,Pi; 
= 2 Res ^^Biq,p,)Biq,p,) 

= 2 Res Res -— ^-^g(P) B{q,pi)B{r,p2) 

g^a r^q [y{q) - y{q)){x{r) - x{q)) 

dE {p\ 

= -2 Res Res -— — B{q,pi)B{r,p2) 

g^a r^q [y[q) - y{q))[x[r) - x[q)) 

= -2 Res ^^^B{q,p,)B{q,p2) 

q^a Uj{q) 

j^gg B{q,p)B{q,pi)B{q,p2) 
q^a dx{q)dy{q) 

where we have used Eq. flA-2p . Eq. flA-ll) and Eq. flA-31) . □ 
Theorem 14.21 proof: 

It is obvious from the definition that if p is away from branchpoints, the residues 
are finite integrals, and wj:^-^ is finite. The only poles can be obtained when p pinches 
an integration contour, i.e. at branch points. 

Then, it is easy to see by recursion on k and g that it holds for any p,. 

□ 



Theorem 14.31 proof: 

The first one is a property of dEq{p), and the second one follows from recursion on 
k and g, and it holds for 1^2°-' = B. n 
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Theorem 14.41 proof: 

The case k = 1, g = comes from Eq. flA-ll) . and by integration: 

Y,dE,{p')=0 (A-5) 

i 

which proves equation l4-23[ Then, it is clear from the iterative definition of wjf^^ for 
k > 1, that the dependance oiWl^^^{pi,p,p2, ■ ■ ■ ,Pk) in p is a sum of integrals involving 
only the following quantities: 

J2 Res ^^^Biq,p) q,,..., q,) (A-6) 

for some gi, . . . , Using equation lA-ll we have to compute: 
^ ^ ^M^B{q,p^) <^l(g, q,,..., q,) 



« 3 



g^aj 00 j (q) 
dEj^gjqi) dx{q)dx{p) (g) 



> Kes T^TT^; W' ^2, • • • , g^j 

^q^a, ujj{q) {x{q)~x{p)y 

^ p^., ^-^i.gfa) dx{q)dx{p) (.,,{9), n^^w'^^^lnn 

2^9^% ujj[q) [x[q) ~ x[p)Y \ 



1 V- V-r>p^ dEU^i) dx{q)dx{p) .,,{9) , j . 
' (A-7) 



where the second equality holds due to IA-21 the third equality holds due to equation 
l4-23[ and the last equality holds because that last expression has no poles at the 
branchpoints. This proves equation 14-241 □ 



Theorem 14.51 proof: 

It is clearly a rational function of x{p) because it is a symmetric sum on all sheets. 
From theorem 14.21 the RHS may have poles at branch points, and/or at the poles of 
some y{p^)i and/or when x{p) = x{pi) for some I. Let us prove that the poles at branch 
points actualy cancel. 

Let us denote in this section: 

ul^\q.q\pK) = Y.Y. W^ti^^PjWt;:l{q,PK/j) + wi%'\q,q\pK). (A-8) 

m=0 JdK 

From theorem 14.41 we have: 

Y.U^^\q\q\pK) 

i 

^°Since it may not be clear which branch point dE and uj refers to, one indicates it with an index. 
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= -U]^\q,q,PK)-U\^\q,q,PK)- iUjf\q,q\pK) + Ui'\q,q\pK)) 
(A -9) 

and from theorem 14.41 and theorem \A.2\ only the first two terms have poles at branch 
points, and thus: 



J2 Res fui'\q,q,p^) + ut\q,q,pK 



= 2Wi'^,{p,pK) (A-10) 
where the last equality holds from the definition of wjf^^. Thus we have: 

(A- 11) 

Then we rewrite it in terms of Pjf^ above: 

Wif,{p,PK) 

= -I E ^%T^ dx{qr 

I ^ n^a-j ujj[q) V 

+2j2yiq')dxiq)W^f,iq\pK)) 

i 

-IYI Res (pt'\x{q),pK) dx{qf + 2y(g)rfx(g)iyifi(g,p;,) 

+22/(g)da;(g)W^ifi(g,PAO 



= _1 ^ Res (pi^\x{qUK)dx{qf 
+2iyiq)-yiq))dxiq)w!:'^,iq,PK) 

{A - 12) 

where we have used theorem 14.41 again. That gives 
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_1 J2 Res ^^^Pt\x{q),pK)dx{qY 



2 ^9-% ujj{q) 
(A -13) 

Let us compute that last integral: 

^Res dEMwi%{q,pK) 

= - Res dEM wi'^,{q,pK) + ^Y.i ^(P' / ^iUl^PK) 

= - Res dEMW^'^,{q,pK) + J2du,{p) I wii,{q,pK) 
= - Res dEj,,{p)wi%{q,pK) 

= -W^Up^pk) (A-14) 

where we have deformed the contour of integration using Riemann bilinear identity, 
and then we have used theorem 14.31 Thus we have: 

y Res ^%Mp(f)(a;(g),p^)rfx(g)2 = 0. (A-15) 
Since this holds for any p, we can write for any m > 0: 

= Res {y{p) - y{a.)) - x(a.))™ Y. '^^^ P^'\x{q),pK) dx^qf 
_ Res Res {y{p) - y{a,)) {x{p) - x(a,))- Pi'\x{q),pK) dx{qf 

q^ai p^q,q Uj[q) 

- Res dx{q) {x{q) - x{a,)r Pt'\x{q),pK) 

{A - 16) 

which proves that P^^\x{q),pK) can have no pole at g = Oj. 
□ 



Theorem 14.61 proof: 

Assume this is proved for h < g, and at g, it is proved for / < k. Then we have: 

Wi'\p^,P2) 

^ ' m=0 
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^—^ q^ai uj[q) 

+ Res Res ^ U g wt'\,' , ,)wt-'-'\-^' 

+2 Wt^^ (g) lyf (g', g', g) + 2 X] l^f (g', g, g) W^^'-'^ (g') 

m=l m'=0 
9 m 

+4E E w^^'^(g^g)v^^^"'^(gO<"'"Hg) + w^?"'H?,9^^,^' 



(A - 17) 



m=l m'=0 



and W2^\p2,pi) is given by the same integral except that the order for computing 
residues is reversed, the residue in g is computed before g'. The difference is thus 
obtained by pushing the contour of q' through the contour of g, and is obtained as the 
residue at g = g'. Notice that only W2^\q, q') has a pole at g = g', all the other W^f'^ 
have no poles at g = g' from theorem 14 .21 Thus we have: 



= E. Res,_, Res,,^,^^5(g,g') 

4 Ei=o Wi""^ mwt""^ (g) + 2Wt'^ {q', q) 
= 2E. Res,_... Res,^q'-^'-^B{q,q')Ul,^\q,q') 
= Res Res ^ ^ 5(g, g') (f/^^^ (g, g') + u!,'^ (g, g': 



(A-18) 



where we have used the notation of theorem I4.5[ 
Similarly, for higher values of k, we find: 

Vf^2+fc (Pi ,P2,Pk)- W^ff^ {p2 ,Pi,Pk) 

- y Res dEq{p^)Biq,p2)-dEqip2)Biq,p,) ^(.) (- 

(A -19) 
Then it gives: 

Wlfl^{pX,P2,PK) - W^2¥fc(P2,Pl,Pi^) 

- y Res d^MB{q.P.)- dE,{p2)B{q^p,) ^(.) 

= yRes'-^dJ'-^ (ut\q,q',p,) + Ut\q,q',PK))) 
(A - 20) 
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and by integrating half of it by parts we get: 
W^i+l (Pi > P2 , Pi^ ) - W!f^^ {P2,Pi,Pk) 

= - Res „ ^jrji,) ^ ^(.) - j ^^_2i) 

+ Res (^)^ 
Now we use the following Lemma: 

Lemma A.l If f{q,q') is localy a bilinear differential near a branchpoint ai, with no 
poles, and symmetric in q and q, and in q' and q' , then: 

Biq,p,)dEM- B{q,p,)dEM ^ ^ 

Proof of the Lemma: 

The residue is a simple pole, and we can use formula lA-31 that gives: 
Res B^'i^P^)dE,{p,)-B{q^p,)dEM^^^^^^ 

B{q,pi)B{q,p2) ~ B{q,p2)B{q,pi) 

= Res . X , . X , . X f{q,q) 

q~*a, uj[q)ay[q)dx[q) 
= (A-23) 

□. 



Using this Lemma, as well as theorem 14.51 we get: 

H^i+i(Pl'P2,Px) - Wi%{p2,Pl,PK) 

= - E. Res,^.. (^(^) _ y^q))dx{q)WtMq,pK) 

+ E. Res ,^a. '^^^(^^)^('^-^^)-^';^^(^^)^('^-^^) Wii, (g, pk) 
= 0. 

□ 



Corollary 14.11 proof: 

For any rational function R{x) with no pole at x{ai) we have: 
Res R{x{p))wi'^M Pi,---,Pk) 

Hi 

= ^ Res R{x{p)) {wtMp^p,, ...,pk) + ...,Pk: 
= " (A-25) 

due to theorem I4.4[ 

For m = 0, 1 compute 

I * n — ^rv 



p— >a 
a. 
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h=0 IcK 

1 Res ^f^2yip)dxip)wt^,ip,pK) + 

+ E E K\iMp^^^tiii,MPK/i) + wI:1-'\p,p,pk) 

h=0 ICK 



g 



+ E E K^p^p^^^tiiUp'P^/^^ + wt2\p,p,PK) 

h=0 IcK 

1^ f x{pjrwt\pK) \ 



+i E ilf. ^ ( - Mp)dx{p)wtUp,PK)+ 

+E E + wt,'\p,p,PK) 

h=0 ICK 

+\ E ^ ( - Mp)dx{p)wtMp,pK)+ 

h=0 ICK 



1 ' f xip,r^w^'\pK) \ 



{A - 26) 
due to theorem 14. 5[ □ 

Theorem 14.71 proof: 

The case k = 1, g = is easy: 



Res $(p)S(pi,p2) = d^p2) = yiPi)dx{pi). (A-27) 
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We prove the theorem by recursion on g and k. Suppose it is proved for all k' for 
g' < ^-1, and for fc' < k-l ii g' = g. We write /sT = {1,. . and K' = {1, . ..,k-l}. 
Then we have from Eq. fl4-9p : 

Res 

= Res Res ( J] W^ti(g,Pj,P.)W^£r^(g,P^7^) 



m=0 JC-ff' 

{A - 28) 

Then we exchange the contours of integration: 



Thus 



Res Res = Res Res + Res Res_ . (A-29) 

Pk^a q^a q^a pt^a q^a Pk^q,q 



Res 

Pk^a 

9 



Res Res ^p,) ( E ^7i(?,P.,Pfc)W^iT^(g,Pxv.) 

q^a Pk^a ^[Q) ^ — — 



m=0 JC-ft:' 
9 



m=0 JcK' 



+ Res Res ( ^ W^:i{q,pj,p,)W^-^-\q,p^,/j) 

"j^a Pk^i^q uj{q) V — ^ — ^ 



m=0 JCX' 

g 



(A - 30) 



m=0 JcK' 



The first term is computed from the recursion hypothesis, and the second term can 
exist only if the correlation function containing pk has poles at pk = q or pk = q, and 
from theorem 14.21 this can happen only if the correlation function containing pk is a 
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Bergmann kernel. That gives: 
ReSp,^a$bfc)M/fc^fi(p,Pi, • • ■,Pk) 

= Res,_ ^ ( EL=o EjcK'i^rn + (j + 1) - 2)W^^f,iq,pj)Wt--\q,pj„/j) 
+ TL=o EjcK'i^ig -m) + ik- j) - 2)W^-liq,pj)wt-"'\q,pK'/j) 
+ {2{g - 1) + A: + 1 - 2)wl^l-'\q,q,pK')) 

+ Res,^. Res,,^,,,$(p,) ^ ( + Wi'\q,pK')Biq,pk)) 

= i2g + k-3) Res,^.^ ( E^_o 

+ Res,_^i^ (2/(g)iyif,(g,pK0 + 
= {2g + k-3)Wk{p,pi...,Pk-i) 

+ WTW) {vi^WliMPK') + wi'\q,pK')) + {y{q) - y(g))W^i^^(g,P;,o) 

= {2g + k - ?>)Wk{p,Pi . . . ,Pk-i) 

= (2(7 + A; - 3) . . . ,Pk-i) - Res, ^^dE,{p) W^^'\q,pK') 

= {2g + k-3) Wk{p,pi . . . ,pk-i) + \ Res,^arf5',,,(p) Wl'\q,pK') 

-\ Res,^^dS-Up)wl;'\q,pK') 
= {2g + k-3) Wk{p,Pi . . . + \ ReSg_arfS',,„(p) Wl^\q,pK') 

- 1 ReSq^^dSq^oip) Wlf'^ {q, pk' ) 
= (2^ + A; - 3) Wk{p,pi . . . ,Pk-i) + Res,^aC?V(p) W^i^^l^-^') 
= (2(? + A; - 3) Wk{p,pi . . . - Res, ^pdS,,o{p) wl:'\q,pK') 

= (2(7 + A; - 3) Wk{p,pi . . . - Res g^pdS^M wl:'\q,pK') 

-du\p) ((1 + kt) §^Wi'\q,PK') -^§^Wi'\q,PK')) 

= {2g + k-3) Wk{p,pi . . . ,Pk-i) - Res g^pdSg,o{p) wl:'\q,pK') 

-du\p) §^wl'\q,pK') 
= {2g + k-3) Wk{p,pi . . . ,pk-i) - Res,^pdS,,oip) W^'\q,pK') 
= {2g + k-2)Wk{p,Pi...,Pk-i) 

(A-31) 

□ 



B Variation of the curve 

In this appendix, we proove the theorems stated in sections ([5]) and section ( 16. ip . 
Lemma 15.11 proof: 
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+2 Res Res .^a. ^ l^(r) ^ /(g, g) 

= E, Res„^ /^n(/(g,g)).(,)-2E, Res,_^. g^fi(g) /(g,g) 

+2 Res Res .^a. ^ ^^(r) ^ /(g, g) 

= Res„ ^ I^n(/(g,g)).(,) Res,^„, Res_,^^%|g^n(r) /(g,g) 

+2 E, Res E. Res .^a. ^ ^^(r) ^ /(g, g) 

= E,Res„^i)n(/(g,g)L(,) 

+2 E. Res r^a. E, Res ^ ^^(r) ^ /(g, g) 

(B-1) 



Theorem 15.11 proof: 

This theorem straightforwardly comes from the diagrammatic rules described in the 
preceding paragraph except for the variation of F^^\ 

Let us prove it for F^^^ with = B(p, q)A{q). One has 



Res, 



Ic dEg{p)A{p)dz,(q)dzi(q) 



where Zi is a local variable near the branch point and Sb is the corresponding 
Bergmann projective connection. Since the last term has a simple pole at the branch 
point ttj, one can write 

_ j^gg jcdEg{p)A{p)dzi{q)dzi{q) l^^/^x ^ -- V Res ^'^^'^^ 

g^a ci;(g) 6 2 dy{ai) q^a,i dx{q) 

= —SnlriTBx- 
(B-S) 

On the other hand, one can express the first term thanks to the local variables Zi 
and compute 

J^dEg{p)A{p)dzi{q)dZi{q) ^^(u 

(5-4) 

provided that Zi{q) — Zi{q) = 2zi{q). □ 



Theorem ( 15.31) proof: 
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Using theorem 14. 7[ we have: 

(2 - 2g)F^3) = _ j^es (^y^ia) j^g.g) 

a 

and using eq J3-63l we can choose for any arbitrary o': 

which imphes: 

(2 - 2g)F^9) = _ Res^^l^^^^ 

= Resp^a Res, ^^V^{q)dSp,o'{q)wi'\p) - Y.Jc. Res ^.^^ j;^^JSp^o'{q)w'f\p) 
- Q Res,^a §,^^^ dSpA<lW['\v). 

(B-7) 

Since the poles a and the branch points Oj do not coincide, one can exchange the order 
of integration. Then, one can move the integration contours for p in order to integrate 
only around the last pole p —>■ q: 

(2 - 2g)F(3) 

= Eatajg^, ReSp^,dSp,o'iq)wi'\p) - E„ Res,^„ ReSp^,V^iq)dSp,o'iq)wi'\p) 
+ Res p^,dSpA<lWi'\p) 



(B-8) 



□ 



Theorem (16.11) proof: 

We have 

2in-^W2°\pi,P2) = -{2inf{dui{pi)duj{p2)+dui{p2)duj{pi)) 

2 



(5-9) 



{B{pi,r)B{p2, s) + B{p2, r)B{pi, s)). 



Thus the theorem holds for k = 2 and g = 0. And by integration we get 
d 

- 2m^—dEq{p) = 2{mf{dUi{p){uj{q) - Uj{q)) + duj{p){ui{q) - Ui{q))) 



2 

[B - 10) 



\l I iBip,r)dE,is)+B{p,s)dE,{r)). 
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IX X R i?(p,r)dE,(s)+i?(p,s)dg,(r) ^ A 

+ Y.w^%{q.pj)wt'Mq,PK/j) 



= I ieB, {BiP^ r)wl'M-',pK) + Bip, .)iyif,(r,p^)^ 

+2W^^i(g, g, pj, r)iyi^_-i7'^ (Pi^/J, ^) + 2W^jS(9, PJ, r)wtjlt\q, Pk/j, s)) 

+ 1 ieBjseB, Ej W^i-iS(?,P/^/j)x 

X (<3 ^, ^) + 2 El r)l^i':-!;Vi(Pi^/(^uL), .)) 

(B-11) 

We regroup together all terms with 2 W^s and all terms with 3 W^s: 
= \ LBjseB, {B{p,r)wtU^,pK) + B{p,s)Wi%{r,pK)) 

+ 1 ^es W^i+f ^ (g, Pk, r, s) 



+ ieB, Lb. f^es 



X 



X ( E. El <1(^,Pl, rWt-'I'UpKKJ^L). s)W^fcS(g,PA7j) 
X ( E. El <1(?,Pl, r)wt^UPKm, s)wt;iMPK/j) 

(B-12) 

We recognize the recursion relation eq l4-9l in lines 2 to 5, and in lines 6 to 8, and this 
gives the theorem. 

The theorem for the free energies, is easy obtained using theorem 14.71 □ 
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C Proof of the symplectic invariance of and F^^) 
C.l fw 

Let us study how F^^^ is changed under the exchange of the roles of x and y. For this 
purpose, we define the images of F*^^) and Wi^\p) under this transormation, i.e. 

and 

FW = -^-lnTBy-^lnl[x'{b,) (C-2) 

j 

where b denotes the set of y-branch points, q is the only point satisfying y{q) = y{q) 
and approaching a branch point bj when q —>■ bj, TBy is the Bergmann tau function 
associated to y and x'{bj) = According to theorem 15. 1[ for any variation Q of 

the curve, the variation of the free energy reads 

SnF^'^ = l^wi'\p)A{p). (C-3) 

Thus, the variation of the difference between the two "free energies" reads 

Sn - = ^ {wi'\p) + Wi'\p)) A(p). (C-4) 

In order to evaluate this quantity, one needs the following lemma: 
Lemma C.l For any choice of variable z: 



Wi'\p) +Wi'\p) = ^d. 



1 / , , X'Y X"2 X'" /// 

X y \ X y X y^ y 



(C-5) 



where the derivatives are taken with respect to the variable z and Sbz denotes the 
Bergmann projective connection associated to z. 

proof: 

From the definitions , one easily derives that 

W['\p) + W['\p) = Res -— p^j^ (C-6) 

{x{p) - x[q))[y[p) - y[q)) 

Let us now expand the integrand in terms of an arbitrary local variable z when 
p q. The different factors read 

- ' +is«.(p)-iM^si^(,)+o((.(rt-.M)') (c-7) 



dz{p)dz{q) {z{p) — z{q)Y 6 12 
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and 



y{p) - viq) (zip) - z{q))y'{p) 

+ {z{p)-z{q)f 



y" 

1 - {z{p) - z{q)) — 

y'" y'" 



^ {z{p)-z{q)f / y- ^ ^y"'y" 3y"' 



24 V y' y'^ y''^ 

+0{{z{p) - z{q)Y) 



Inserting it altogether inside Eq. (lC-61) . one can explicitly compute the residues and 
one recognizes the formula IC-5[ □ 

We can now prove the following theorem stating the symmetry of F^^^ under the 
exchange of x ^ y: 

Lemma C.2 The two free energies transform in the same way under any variation of 
the moduli of the curve: 

6nF^'^ = 6nF^'^ (C-9) 

proof: 

We already know that 

(C-IO) 

for an arbitrary variable z. 

We just have to check that this quantity vanishes for the transformations corre- 
sponding to varying the moduli of the curve. Because the function inside the differ- 
ential wrt p vanishes at the poles of ydx, one can check that it is indeed the case. 
□ 

Thus the first correction to the free energy satisfies the following variation under 
symplectic transformations: 

Theorem C.l F'^^'^ does not change under the following transformations of 8: 

• y y + R{x) where R is a rational function. 

• y ^ cy and x ^ ^x where c is a non-zero complex number. 

• X ^ X + R{y) where R is a rational function. 

• y ^ X and x y. 

F^^^ is shifted by a multiple of in/ 12 when £ is changed by x ^ —x. 



80 



proof: 

The first transformation is obvious from the definition since neither Ititex nor y'{ai) 
changes. 

The second one follows from theorem 2 in p7] for f = x and g = ^ which shows 
that the variations of ln{TBx) and y'{ai) compensate. 

The fourth one is nothing but lemma (10.20 and it gives the third one when combined 
with the first. 

The last transformation holds because Iutbx is left unchanged and y'{ai) changes 
sign. 

□ 

C.2 

Theorem C.2 F*^'') does not change under the following transformations of 8: 

• y y + V{x) where V is a polynomial. 

• y —>■ y and x —x. 

• y ^ cy and x ^ ■^x where c is a non-zero complex number. 

• X X + V{y) where V is a polynomial. 

• y ^ X and x y. 

proof: 

The second and third transformations come from the structure of F^^^ which is 
bilinear in objects proportionnal to ydx. 

One can obtain the last one following the same method as for F^^^: we compute 
the variation of the difference of the two free energies under the changes of the moduli 
of S and show that they vanish. 

Let us now show the first invariance. In this case, the variation of the free energy 
is S'p(x)dxF^^\ Since V{x) is a polynomial one can write it 

V{x{p)) = Res B{p,q)Q{x{q)) = - V Res B{p,q)Q{x{q)), (0-11) 

a 

where Q{x) is also a polynomial in x. Then, theorem (15. ID implies the invariance of 
F^^^ under this transformation. 

The fourth transformation is a combination of the first and the last one. 

□ 
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D Matrix model with an external field 



We consider the matrix model with external field defined in section fll0.4l) : 

Z{K) := [ rfAfe-^^^(^(^^)-*^^) (D-1) 




J Hn 

where we assume that A is the diagonal matrix: 

A = diag ( Ai, . . . , Ai, A2, . . . , A2, . . . , A^, . . . , As ) (D-2) 

and V'{x) is a rational fraction with denominator D{x): V'{x) = ^^^^^^j—- 

In particular, the polynomial S{y) := 11^=1 (?/ ^ ■^i) minimal polynomial of A. 

We define the correlation functions uJfc(xi, . . . ,Xk) := N^~'^ (Y[i=i x ~ai ^ ^"^^ 

their l/N"^ expansion 

Wk{^K) = (D-3) 

h=0 

We also define the auxiliary functions 



and 



Notice that UkX^j Vi ^k) is a polynomial in y of degree s — 1, and D{x)Pk{x, y; xj^) 
is a polynomial in x of degree d—1 and in y of degree s — 1 (note that Pq corresponds 
to P in Eq. ([lESHD). 

It is convenient to renormalize those functions, and define: 

Uk{x, y; xk) := Uk{x, y; xr) - Sk,oS{y) (D-6) 

and ^ 

Wfc(xi^) := MJfc(xii-) + ^^^-T^. (D-7) 

(xi - X2y 

D.l Loop equations 

Consider the change of variables 

SM= (D-8) 

x-M y-A ^ ' 
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You get the loop equation 

Wi{x)uo{x,y) + ^Ui{x,y]x) = V'{x)uo{x,y)-PQ{x,y)-yuo{x,y) + S{y)wi{x) (D-9) 
i.e. 

{y+wi{x)-V\x)){Mx,y)-S{y)) + ^Mx,y-,x) = {V'{x)-y)S{y)-Po{x,y). (D-10) 
We define the polynomial both in x and y 

EmsA^, y) ■■= iiV\x) - y)S{y) - Po{x, y)) D{x) (D-11) 

and 

Y{x):^V'{x)-wi{x). (D-12) 
The loop equation thus implies: 

1 

{y - Y{x))uo{x, y)D{x) + -^-"1(3;, y; x)D{x) = EMext{x, y) (D-13) 
and in particular 

EMeAx.Y{x)) = ■^mix,Y{x);x)D{x). (D-14) 
The leading order of the topological expansion reads 

E^SeJ^,Y(x))^0 (D-15) 

which defines an algebraic curve. 

D.2 Leading order algebraic curve 

Let us study the curve SMext{x,y) — E^^^^^(x,y) — defining a compact Riemann 
surface E and two functions x and y defined on it. 

Because y is a solution of a degree s + 1 equation, SMexti^, v) s + 1 x-sheets. 
The sheets can be identified by their large x behavior: 

• in the physical sheet, we have Y{x) ~ V\x) — 1/x + 0{l/x'^) 

• in the other sheets, Y{x) + + 

Let us note by e E with i — . . . s the different points of the curve whose 
x-projection are i.e. 

Vi,i x{p')^x{p'). (D-16) 

The superscript corresponds to the point in the physical sheet. 

Prom the correlation functions previously defined on the x and y projections, one 
defines the corresponding meromorphic 1-forms on the curve as follows: 

W^fc(PK) := w^fc(x(pK)) dx{pi) . . . dx{pk) (D-17) 
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and 

Uk{p, y; Pk) := Uk{x{p),y; x(pk)) dx{p)dx{pi) . . . dx{pk) (D-18) 
as well as there topological expansions 

oo oo 

W^.(pK) = 5^Ar2-2'^W^f(pj,) and f/,(p, Pk) = J] iV^"^'^ f/f (p, 2/; Pk)- 

h=Q h=0 

(D-19) 

D.2.1 Filling fractions and genus 

The curve has a genus g < ds — 1 and we work with fixed filling fractions 



€j := ^ (p ydx. (D-20) 
'Ai 



2m 



D.2.2 Subleading loop equations 

Consider the topological expansion of the loop equation Eq. (]D-14|) . It reads, for /i > 1: 



E^^\x,y) = D{x){y - Y{x))uf\x,y) + D{x)w'^l'),{x)u^^\x,y) 

+D{x) ^w'y^\x)u''Q'''^\x,y) + D{x)uf'^\x,y;x), 

m=l 

{D-21) 

where E^^\x^y) is the /I'th term in the /i^-expansion of the spectral curve. 

D.3 Diagrammatic rules for the correlation functions and the 
free energy 

In this section, one proves that the correlation functions' and the free energy's topologi- 
cal expansion of this model do coincide with the W^lf'^s and F'-'^-^'s defined following the 
definitions of Eq. (14-151) and Eq. (I4-I6p for the classical spectral curve SMext{x,y) = 0. 

D.3.1 The semi-classical spectral curve 

Let us reexpress the semi-classical spectral curve (i.e. the whole formal series 
EMext{x,y)) in terms of the classical one E^^^^^{x,y). 

Theorem D.l 

EMext{x,y) = -D{x) " (^n!=o(l/ - y'{x{p)) + F ^.(pwVa/ )) (D-22) 
= D{x)[{V'{x)-y)S{y)-P^{x,y)] 
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and 



Uo{p,y) = -" (fl^y - V'ixip)) + Ixr ^py}^)^ (D-23) 



where " < . >" means that one replace W2 by W2 in the expansion. 
proof: 

One proves that the j^-expansions of 



E{x, y) = -D{x) (\[{y - V\x{p)) + 1 Tr ^py]^) ) (D-24) 



,i=0 

and 



f/(p, y) = -D{x) ( Hiy - V\x{p)) + 1 Tr ^py^) ) (D-25) 



I 4 = 1 



coincide with the expansion of EMext{x,y) and Uo{p,y). 
Let the topological expansions be 

Eix,y) = J2N-''E^'\^,y) , U{p,y) = Y,N-^m^>^\p,y). (D-26) 



Expanding the expressions of E{x,y) and U{p,y) into cumulants, one recovers 

= {y- Y{x))D{x)ut\x,y) + D{x)w^^\x)U^^\x,y) 

+D{x) w^^\x)ut"'\x, y) + D{x)Uf'~'\x, y- x), 

m=l 

{D - 27) 



which coincides with Eq. (lD-2ip . 

One easily proves that this system of equations admits a unique solution thanks to 
the polynomial properties of U{p, y) and that the leading orders h = coincide. The 
proof is extremely similar to that for the 2-matrix model (cf. theorem 1 in [H]). □ 

D.3.2 Diagrammatic solution 

One has: 

wi'\pi,P2) = B{p,,p2) (D-28) 

where 5 is the Bergmann kernel of the algebraic curve Smsxi- 
The coefficient of of Eq. flD-22p . divided by D{x), is: 



V'{x)+J2^^ = Y.^iP')- (D-29) 
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It implies that 

dx{p)dx{q) 



{x{p) -x{q)y 



J2 Wi'\p\ <D and V/i > 1 , 5^ W^''\p\ q) = 0, (D-30) 



I.e. 



^2(p,g) + 5^^2(p\g) = 0. (D-31) 

i=l 

The coefficient of y"^'^ is: 

(D-32) 



j<j j i<j ^ ' 



Notice that: 



2 

{D - 33) 



i<j ^ ^ 

^-{V'{x) + E A.)' - ^ E (^(P^)' + . 



Thus: 



i ^ ' i ^ ^ 

Notice that the LHS is the ratio of a polynomial in x and D{x): §|| = V'{xf + 

\- A2_ 2 / V'(x)-y'(Afi) \ 

The topological expansion of this equation reads for /i > 1 



From now on, following the lines of [19], one multiplies these equations by | y{J)-^^) ■ 



(D-35) 

(P)- 

takes the residues when p —>■ fia and sums over all the branch points and obtains: 

^1,0 (9) = E R;^!^ r.r.^-.r^^w.r.^ • ^^'^e) 
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Differentiating wrt the potential V{xi), one can finally write down an expression 
for the correlation functions: 



W, 



(h) 



Res, 



{y(p)-y(p)) dx{p) 



(D-37) 



This coincides with the reccursive definition 14-151 and ensures the equality of the 
correlation functions with the former defined special "loop functions". 

Keeping on following p^, one finds that the topological expansion of the free energy 
also coincides with the special free energies defined on 14- 16^ that is the r-function of 
the algebraic curve. 
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